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Abstract 


A search for a high-mass Higgs boson H is performed in the H WW Mv and 
H —> WW —> (vqq decay channels using pp collision data corresponding to an integrated 
luminosity of 20.3 fb”' collected at = 8 TeV by the ATLAS detector at the Large Hadron 
Collider. No evidence of a high-mass Higgs boson is found. Limits on cth x BR(//^WW) 
as a function of the Higgs boson mass mn are determined in three different scenarios: one 
in which the heavy Higgs boson has a narrow width compared to the experimental resol¬ 
ution, one for a width increasing with the boson mass and modeled by the complex-pole 
scheme following the same behavior as in the Standard Model, and one for intermediate 
widths. The upper range of the search is uih = 1500 GeV for the narrow-width scenario 
and iuh - 1000 GeV for the other two scenarios. The lower edge of the search range is 
200-300 GeV and depends on the analysis channel and search scenario. For each signal in¬ 
terpretation, individual and combined limits from the two WW decay channels are presented. 
At niH = 1500 GeV, the highest-mass point tested, cth x BR(// ^ WW) for a narrow-width 
Higgs boson is constrained to be less than 22 fb and 6.6 fb at 95% CL for the gluon fusion 
and vector-boson fusion production modes, respectively. 
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1 Introduction 


The boson discovered in 2012 by the ATLAS [1] and CMS [2] collaborations at the LHC matches the 
predictions for a Standard Model (SM) Higgs boson within the precision of current measurements [3, 4]. 
Several extensions of the SM predict heavy neutral scalars in addition to a low-mass scalar compatible with 
the discovered boson. Examples include generic models in which a low-mass Higgs boson mixes with a 
heavy electroweak singlet [5-10] to complete the unitarisation of WW scattering at high energies. 

This paper reports the results of a search for a heavy neutral scalar by the ATLAS Collaboration in the 
decay mode into two W bosons. Two final states are used: H —> WW —> Mv and H WW —> ivqq {{ = 
e,ii). In these final slates, ATLAS has previously reported Ihe resulls of searches for heavy Higgs bosons 
using 4.7 fb“^ of dala collected al a cenlre-of-mass energy of 7 TeV [11, 12]. In Ihe H WW —> Mv final 
slate, a SM Higgs boson in Ihe mass range 133 GeV < mh < 261 GeV was excluded af 95% confidence 
level (CL), while Ihe H WW —> £vqq final slafe was nol sensitive lo a SM Higgs boson of any mass 
wilh Ihe ^fs -1 TeV dalasel. The CMS Collaboralion has performed a search for a heavy Higgs boson in 
Ihe H^WW and H^ZZ channels [13]. From a combination of Ihe Iwo channels, a hypolhelical second 
Higgs boson wifh couplings idenfical lo fhose predicled by fhe Slandard Model is excluded in Ihe mass 
range 145 GeV < mn < 1000 GeV. 

The analyses reported here improve Ihe resulls in Refs. [1 1, 12] by using an integrated luminosity corres¬ 
ponding lo 20.3 fb“^ of pp collision dala al = 8 TeV collected by Ihe ATLAS detector. Bolh analyses 
are designed to be sensitive to a heavy Higgs boson produced Ihrough eilher or bolh of Ihe gluon-fusion 
(ggF) or vector-boson fusion (VBF) processes. Bolh also use a profile-likelihood til to a dislribulion in 
which Ihe hypolhelical signal is peaked bul background is monolonically decreasing in Ihe search range in 
order to lesl for Ihe presence of signal. The H WW —> Mv analysis uses Ihe dileplon Iransverse mass 
dislribulion for Ihe discriminanl because Ihe Iwo neulrinos in Ihe final slate resull in insufficienl kinematic 
information to reconslrucl Ihe invarianl mass of Ihe WW system. The H —> WW —> ivqq analysis uses 
as Ihe discriminanl Ihe invarianl mass of Ihe WW system, reconslrucled using Ihe W mass as a kinematic 
conslrainl to recover Ihe neulrino momenlum up to a Iwofold ambiguity. The resulls of Ihe searches are 
interpreted in Ihree scenarios: 

1. A Higgs boson wilh Ihe couplings predicted by Ihe SM for a Higgs boson al high mass and a 
widlh correspondingly increasing wilh mh, and Ihe lineshape modeled by Ihe complex-pole scheme 
(CPS) for mosl mass hypolheses, as explained in Sec. 2. Accordingly, Ihis is referred to as Ihe CPS 
scenario. 

2. A Higgs boson wilh a narrow widlh: labelled as narrow-widlh approximation (‘NWA’). 

3. An intermediate-widlh (TW’) scenario, motivated by Ihe eleclroweak singlel model. 

Section 2 of Ihis paper discusses Ihe CPS lineshape model. Section 3 describes Ihe ATFAS detector, 
Ihe dala sample and physics objecl reconslruclion. Section 4 summarises Ihe simulation of signal and 
background samples. The evenl selection and background estimation techniques used in Ihe analyses are 
described in Sections 5 and 6. Systematic uncerlainlies affecting Ihe analyses are discussed in Section 7. 
Dislribulions of Ihe discriminanls are shown in Section 8. Section 9 presenls Ihe inlerprelalions of Ihe 
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results from the H —> WW —> Mv and H —> WW —> ivqq final states, as well as from their combination, 
in the scenarios listed above. Conclusions of the study are given in Section 10. 


2 CPS lineshape model for a heavy Higgs boson 

Narrower widths are allowed in general for Higgs bosons in extensions to the Standard Model, but to 
explore the implications of the width of the additional Higgs boson, the data are also interpreted using a 
signal hypothesis with a lineshape and width identical to a SM Higgs boson. The width of a SM Higgs 
boson increases with increasing mass. For example, it is ~ 30GeV at mn - 400 GeV, and increases to 
~ 650 GeV at mn - 1000 GeV. Up to mn ~ 400 GeV, the lineshape of the WW invariant mass (mww) 
distribution is well described by a Breit-Wigner distribution with a running width, meaning that the Higgs 
boson propagator is calculated for each event based on mww as described in Ref. [14]. For nin > 400 GeV, 
the complex-pole scheme [15-17] provides a more accurate description. The CPS propagator is therefore 
used to describe the lineshape of the Higgs boson produced via both the ggF and VBF processes for 
niH > 400 GeV [18-20]. The limits using this signal hypothesis are labeled “CPS scenario” even though a 
Breit-Wigner distribution is used for m// < 400 GeV. For that mass range the distributions are similar, so 
this is a minor simplification. 

For a Higgs boson with a large width, the production cross section as well as the shapes of kinematic vari¬ 
ables are affected by the interference between signal and non-resonant WW background. The interference 
is small for itih < 400 GeV, but is significant at higher masses, since it increases with increasing Higgs 
boson width. The effect of the interference is included in the signal samples which use the CPS lineshape, 
i.e. niH > 400 GeV. The interference calculations are described in Section 4. 


3 Data sample and object reconstruction 


The ATLAS detector [2 1 ] is a general-purpose particle detector used to investigate a broad range of physics 
processes. It includes inner tracking devices surrounded by a superconducting solenoid, electromagnetic 
(EM) and hadronic calorimeters and a muon spectrometer with a toroidal magnetic field. The inner de- 
tecfor (ID) consists of a silicon pixel detector, a silicon microstrip detector, and a straw tube tracker that 
also has transition radiation detection capability. The ID provides precision tracking of charged particles 
with pseudorapidity ^ \q\ < 2.5. The calorimeter system covers the pseudorapidity range \q\ < 4.9. It is 
composed of sampling calorimeters with either liquid argon or scintillator tiles as the active medium. The 
muon spectrometer provides muon identification and measurement for \q\ < 2.7. During Run 1 of the 
LHC, the ATLAS detector used a three-level trigger system to select events for offline analysis. 

Owing to the high LHC luminosity and a bunch separation of 50 ns, the number of proton-proton interac¬ 
tions occurring in the same bunch crossing is large (on average 20.7 in 2012). Proton-proton interactions 

' ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the i/-axis points 
upward. Cylindrical coordinates (r, (j)) are used in the transverse plane, being the azimuthal angle around the beam pipe. The 
pseudorapidity is defined in terms of the polar angle 6 as 77 = - lntan(0/2). 
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Table 1: The minimum transverse momentum (pj) requirements used at the different levels of the trigger. An “i” 
next to the threshold value indicates an isolation requirement that is less restrictive than the isolation requirement 
used in the offline selection. The single-lepton triggers with higher-px thresholds are more efficient at high lepton 
Pt than the lower-px triggers because of this isolation requirement. For dilepton triggers, the pair of thresholds 
corresponds to the leading and subleading lepton, respectively. The 0 GeV in the line describing the dimuon trigger 
indicates that only one muon is required at Level 1. 


Name 

Level-1 trigger 

High-level (software) trigger 

Single lepton 

e 

18 OR 30 GeV 

24i OR 60 GeV 


15 GeV 

24i OR 36 GeV 

Dilepton 

e, e 

10 AND 10 GeV 

12 AND 12 GeV 

p,p 

15 ANDO GeV 

18 AND 8 GeV 

e,p 

10 AND 6 GeV 

12 AND 8 GeV 


in nearby bunch crossings also affect the detector response. These additional interactions are collectively 
referred to as event “pile-up”^ and require the use of dedicated algorithms and corrections to mitigate its 
effect on particle identification, energy calibrations, and event reconstruction. 

The triggers used in these analyses are listed in Table 1, together with the minimum transverse momentum 
(Pt) requirements at the different levels. Both the H —> WW —> Mv and the H WW —> ivqq analyses 
use the single-lepton triggers while the dilepton triggers are used only by the H WW £v£v analysis. 
The lepton trigger efficiencies are measured using Z boson candidates as a function of lepton pj and rj. 
The single-lepton trigger efficiencies are approximately 70% for muons with | p | < 1.05, 90% for muons in 
the range 1.05 < | p | < 2.40, and >95% for electrons in the range 1 77 1 < 2.40. Dilepton triggers increase the 
signal acceptance for the H —> WW —> Mv analysis by enabling lower lepton px thresholds to be used. 

Events are required to have a primary vertex consistent with the known interaction region, with at least 
three associated tracks with px > 0.4 GeV. If multiple collision vertices are reconstructed, the vertex with 
the largest summed p^ of the associated tracks is selected as the primary vertex. Data quality criteria 
are applied to events to suppress non-collision backgrounds such as cosmic-ray muons, beam-related 
backgrounds or noise in the calorimeters. The resulting integrated luminosity is 20.3 fb“^ at = 8 
TeV. 

Electron candidates are required to have a well-reconstructed track in the ID pointing to a cluster of 
cells with energy depositions in the EM calorimeter. They are required to be in the range | p | < 2.47, 
excluding the range 1.37<|p|<1.52 which corresponds to the transition region between the barrel and 
the endcap calorimeters. Only electrons with Ej > 15 GeV are used in the analysis. The fine lateral 
and longitudinal segmentation of the calorimeter and the transition radiation detection capability of the 
ID allow for robust electron reconstruction and identification in the high pile-up environment. Criteria 
including the calorimeter shower shape, the quality of the match between the track and the cluster, and 
the amount of transition radiation emitted in the ID, are used to define a set of identification criteria [22- 

^ Multiple pp collisions occurring in the same (nearby) bunch crossing are referred to as in-time (out-of-time) pile-up. 
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24]. The “tight” criteria, which have the best background rejection, are used in the H WW —> ivqq 
analysis. The H WW —> Mv analysis uses the “medium” selection, which is more efficient but admits 
more background, for electrons with Et > 25 GeV. For electrons with 15 GeV < Et < 25 GeV, a 
likelihood-based electron selection at the “very tight” operating point is used for its improved background 
rejection. 

Muon candidates are identified by matching tracks reconstructed in the ID with tracks reconstructed in 
the muon spectrometer [25]. The muon spectrometer track is required to have a track segment in each 
of the three layers of the spectrometer, while the ID track must have a minimum number of associated 
hits in each subdetector. In the H —> WW —> Mv analysis, muons are required to have 1 77 ! <2.5 and 
Pt > 15 GeV. For the H WW — > £vqq analysis, muons must satsify \ q\< 2.4 and pj > 25 GeV, since 
the sole lepton in the event must be within the acceptance of the trigger. 

Additional selection criteria on the lepton isolation and impact parameter are used to reduce backgrounds 
from non-prompt leptons and lepton-like signatures produced by hadronic activity. These requirements 
are identical for the H —> WW —> Mv and H WW —> ivqq analyses. Lepton isolation is defined using 
track-based and calorimeter-based quantities. The track isolation is based on the scalar sum Epj of all 
tracks with pj > 0.4 GeV in a cone in q-(p space around the lepton, excluding the lepton track. The 
cone size is AR = + (Aq)^ = 0.3. The track isolation requires that XpT divided by the electron 

transverse energy Ej (muon pj) be less than 0.10 (0.12) for £'t(Ft) > 20 GeV. For electrons (muons) 
with 15 GeV < £'t(Ft) < 20 GeV, the threshold is 0.08. 

The calorimeter isolation selection criterion is also based on a ratio. For electrons, it is computed as the 
sum of the transverse energies, l,Ej, of surrounding energy deposits (topological clusters) in the EM and 
hadronic calorimeters inside a cone of size AR = 0.3 around the candidate electron cluster, divided by the 
electron Ej. The cells within q x (f> = 0.125 x 0.175 around the cluster barycentre are excluded. The 
pile-up and underlying event contribution to the calorimeter isolation is estimated and subtracted event- 
by-event [26]. Electrons with Ej > 20 GeV are required to have relative calorimeter isolation less than 
0.28. Eor 15 GeV < Ej < 20 GeV, the threshold decreases to 0.24. 

Eor muons, the relative calorimeter isolation discriminant is defined as EEj of EM and hadronic calor¬ 
imeter cells above a noise threshold inside a cone of size AR = 0.3 around the muon direction divided 
by the muon pj. All calorimeter cells within a cone of size AR = 0.05 around the muon candidate are 
excluded from the sum. A correction based on the number of reconstructed primary vertices in the event 
is applied to 'LEj to compensate for extra energy due to pile-up. Muons with pj > 25 GeV are required 
to have relative calorimeter isolation less than 0.30. Below that pj value the threshold decreases in steps 
with decreasing pj, with a minimum value of 0 . 12 . 

The significance of the transverse impact parameter, defined as fhe transverse impact parameter do di¬ 
vided by its estimated uncertainty, o-^q, of tracks with respect to the primary vertex is required to sat¬ 
isfy Idol/cr^g < 3.0. The longitudinal impact parameter zo must be |zol sin0 < 0.4 mm for electrons and 
Izol sin 0 < 1.0 mm for muons. 

Jets are reconstructed from topological clusters of calorimeter cells [27-29] using the anti-kf algorithm 
with a radius parameter of 0.4 [30]. The jet energy dependence on pile-up is mitigated by applying two 
data-derived corrections. One is based on the product of the event pj density and the jet area [26]. The 
second correction depends on the number of reconstructed primary vertices and the mean number of 
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expected interactions. After these corrections, an energy- and 77 -dependent calibration is applied to all 
jets. Finally, a residual correction from in situ measurements is applied to refine the jet calibration. In 
both analyses, jets are required to have px > 25 GeV if they have 1 77 1 < 2.4. For jets with 2.4 < 1 77 1 < 4.5, the 
px threshold is raised to 30 GeV. The increased threshold in the forward region reduces the contribution 
from jet candidates produced by pile-up. To reduce the pile-up contribution further, jets within the inner 
detector acceptance are required to have more than 50% of the sum of the scalar px of their associated 
tracks due to tracks coming from the primary vertex. 

Very heavy Higgs bosons give large momenta to their decay products. In the H WW —> €vqq analysis, 
the dijet system produced by the W boson from such a decay is highly boosted and the jets overlap in 
the calorimeter, so they cannot always be resolved with the standard anti-kf algorithm. Therefore, in this 
analysis the hadronic W decay can also be reconstructed as a single jet found by the Cambridge/Aachen 
algorithm [31], built from topological clusters with a radius parameter of 1.2, referred to as large-R jets. 
These jets can mitigate the loss of signal efficiency, and background can be reduced by selecting those 
with features typical of jets originating from two hard partons. These jets are selected using a mass-drop 
filter algorithm [32]. 

Jets containing f 7 -hadrons are identified using a multivariate Z 7 -tagging algorithm [33, 34] which combines 
impact parameter information of tracks and the reconstruction of charm and bottom hadron decays. These 
analyses use a working point with an efficiency of 85% for f 7 -jets and a mis-tag rate for light-flavour jets 
of 10.3% in simulated tt events. High b-]&t tagging efficiency maximises top-quark background rejection, 
which is important for the sensitivity of analysis categories that require one or more jets. 

In the H WW £v£v analysis, two different definitions of missing transverse momentum are used. 
The calorimeter-based definition, is the magnitude of the negative vector sum of the transverse 

momenta of muons, electrons, photons, and jets. Clusters of calibrated calorimeter cells that are not asso¬ 
ciated with any of these objects are also included [35]. This definition takes advantage of the hermeticity 
of the calorimeters and their ability to measure energies of neutral particles. However, the resolution of the 
calorimeter-based quantity is degraded by the significant event pile-up. The resolution can be improved 
by using track-based measurements of the momenta of particles not associated with an identified object to 
replace the calorimeter cell based measurements. The tracks are required to have px > 0.5 GeV and must 
originate from the primary vertex. In practice, the px of these tracks replace the Ej of calorimeter cells 
not associated with identified objecfs. The accurafe primary-vertex association makes the track-based 
measurement more robust against pile-up than the calorimeter-based measurement. The quantity thus 
formulated is referred to as p™®®. 

Using the direction of p™'®** relative to leptons and jets improves the rejection of Drell-Yan backgrounds 
in the H —> WW —> Mv final stale. A quantity p™**] is defined as follows: 



p™®® sinA 0 , 


'near 


if A0near ^ 7r/2 
Otherwise, 


'near 


( 1 ) 


,miss 


,miss 


where A^near is the azimuthal distance of the and the nearest high-px lepton or jet. A calorimeter- 
based quantity is defined similarly. In Drell-Yan events, in which F™*® arises from mismeasurement 
of the Et or px of objects, these quantities tend to have small values, while in events with genuine F™*® 


they have larger values on average. Selection using these quantities therefore rejects Drell-Yan events in 
preference to signal events. 


4 Signal and background simulation 

This section describes the signal and background Monte Carlo (MC) generators used in the analyses, the 
different signal models used in the hypothesis tests, and the cross-section calculations used to normalise 
backgrounds. 

For most processes, separate MC programs are used to generate the hard scattering and to model the parton 
showering (PS), hadronisation, and underlying event (UE). PythiaS [36], Pythia6 [37], Herwig [38] and 
Sherpa [39] are used for the latter three steps for the signal and for some of the background processes. 
When Herwig is used for the hadronisation and PS, the UE is modelled using Jimmy [40]. 

The parton distribution function (PDE) set from CTIO [41] is used for the Powheg [42] and Sherpa 
samples, while CTEQ6E1 [43] is used for the Alpgen [44], Herwig, GG2WW [45], Pythia6 and PythiaS 
samples. Acceptances and efficiencies are obtained from a full simulation [46] of the ATEAS detector 
using either Geant4 [47], or Geant4 combined with a parameterised calorimeter simulation [48]. The 
simulation incorporates a model of the event pile-up conditions in the data, including both in-time and 
out-of-time pile-up. 


4.1 Simulation and normalisation of signal processes 

The Powheg generator combined with PythiaS is used to model all signal processes. Heavy Higgs 
boson production via the ggp and VBE processes are considered in both the H WW —> €v£v and 
H WW —> £vqq analysis channels. Contributions from Higgs-strahlung and ttH production mechan¬ 
isms are not considered owing to their very small cross sections at high Higgs boson masses. Eor leptonic 
W decays, the small contribution from leptonic IT —> rv —> twv decays is included. 

The ggp signal cross-section calculation includes corrections up to next-to-next-to-leading order (NNEO) 
in QCD [49-54]. Next-to-leading-order (NEO) electroweak (EW) corrections are also applied [55, 56], 
as well as QCD soft-gluon resummations up to next-to-next-to-leading logarithmic order (NNEE) [57]. 
These calculations are described in Refs. [14, 58, 59] and assume factorisation between the QCD and EW 
corrections. The VBE signal cross section is computed with approximate NNEO QCD corrections [60] 
and full NEO QCD and EW corrections [61-63]. The total width for the CPS scenario follows the SM 
predictions for high mass and has been calculated using Hdecay [64]. The branching fractions for the 
decay to WW as a function of nin have been calculated using Prophecy4e [65, 66]. 


4.1.1 Signal samples for CPS scenario 

Simulated Higgs boson samples with the width predicted by the SM as a function of mu are generated 
using PowHEG-i- PythiaS, at 20GeV intervals for 220 GeV < mu < 580 GeV, and at 50GeV intervals for 
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600 GeV < niH < 1000 GeV. The CPS-scenario interpretation is not performed for mn > 1000 GeV be¬ 
cause of the large width of the resonance. For mu < 400 GeV, ggF and VBF samples are generated with 
the running-width Breit-Wigner propagator described in Section 2. For mu > 400 GeV, samples are gen¬ 
erated using a CPS propagator. The calculations using the Breit-Wigner and the CPS propagators are in 
good agreement in the mass range below 400 GeV. 

Calculations of the interference effect between resonant and non-resonant gg —> WW production are 
available only at leading-order (LO) accuracy in QCD. Therefore, this effect is not directly included in 
the generation of the ggF and VBF CPS-scenario signal samples, and is implemented via event weighting 
at particle level. The full weighting procedure, including the treatment of associated uncertainties, is 
described in detail in Ref. [10] and summarised here. 

For ggF signal samples, the interference weights are computed at LO using the MCFM [67] program, and 
rescaled to NNLO following the recommendations given in Ref. [10]. EW corrections are also included 
in the NNLO result used in the rescaling. The interference changes the total cross section. For mu > 
400 GeV, it increases with increasing rnn, with an enhancement of almost a factor of four for mn - 
1 TeV [18]. The interference is negative below ~ 400 GeV, but changes the cross section by 10% or 
less. The weighting procedure has also been performed with the GG2WW program; the results show good 
agreement with those using MCFM. The procedure accounts for theoretical uncertainties associated with 
the LO-to-NNLO scaling as well as those due to missing higher-order terms in the presently available 
interference estimation. The weights are applied to the signal samples only, because in the absence of 
signal there is no effect on the background. The sum of the weighted signal and the continuum WW 
background spectra approximately reproduces the results of the full calculation. 

For VBF signal samples, the REPOLO tool provided by the authors of VBENLO [68] is used to extract 
the interference weights. QCD scale and modelling uncertainties associated with the weights are also 
estimated using REPOLO. In this case, the LO-to-N(N)LO differences are expected to be small [60-62, 
69, 70], and no explicit uncertainty is assigned to take these differences into account. Because not all 
of the information needed for the weight calculation is present in the fully reconstructed Monte Carlo 
samples, the weights are parameterised as a function of m^w and m//. A closure test comparing the signal 
lineshapes produced by the reweighting compared to the full calculation for the interference effects shows 
some differences, which are largest for mww far from mn, but do not exceed 10%. These differences are 
treafed as a systematic uncertainty on the signal. 

Eor both ggE and VBE signal, the weights accounting for interference effects are calculated for each 
Higgs boson mass at which the samples are simulated, and applied as a function of mww in the range 0.5 
< mwwI^H < 1-5. The procedure modifies the event kinematics, including the ott distribution used in 
the H WW —> Mv analysis. It has been shown that the weights describe the effect of interference on 
all kinematic variables used in the analyses [10]. 


4.1.2 Narrow-width signal samples 

Eor the narrow-width Higgs boson scenario, signal samples are generated with Powheg-i-Pythia8 using a 
fixed 4.07 MeV-wide Breit-Wigner lineshape at 100 GeV intervals for 300 GeV <mH< 1500 GeV. Owing 
to the small width, the effect of interference between signal and continuum background is negligible over 
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the full mass range explored in the analyses [18, 71], therefore no interference weights are applied to these 
samples. 


4.1.3 Signal samples for intermediate-width scenario 

The intermediate-width scenario signal samples are derived by weighting the CPS signal samples to 
modify the width and lineshape and to account for interference. The lineshape of the heavy Higgs bo¬ 
son is weighted to one derived from a running-width Breit-Wigner propagator, and to scale the width 
down from the SM width. The interference weights are derived using the MCFM and REPOLO tools 
respectively for ggF and VBF signals, as in the CPS scenario, and are computed as a function of the mod¬ 
ified width of the heavy scalar. The interference is a significant effect for T// > lOGeV. The weights are 
applied to the distribution and modify the event kinematics accordingly. 

Intermediate-width signal scenarios are explored for a mass mn between 200 GeV and 1000 GeV and a 
width in the range 0.2r//,sM < T// < 0.8^//_sM^ where T// is the width of the hypothetical particle and 
r H,SM is the width of a SM Higgs boson for the same mass. The extremes of a very narrow width and the 
same width as the SM are covered by the NWA and CPS scenarios. 

4.2 Background processes 

4.2.1 Background processes for the H -» WW -» (v(v analysis 

The MC generators used to simulate the background processes in the H WW —> £v£v analysis, and the 
cross sections used to normalise them, are listed in Table 2. In this table, all W and Z boson decays into 
leptons (e, fj., r) are included in the corresponding products of the cross sections (cr) and the branching 
ratios (BR). 

Cross sections for top-quark and diboson processes are computed as follows. The tt production cross sec¬ 
tion is normalised to the NNLO-i-NNLL computation from TOP-i-i-2.0 [72-74], and single-top processes 
are normalised to NNLL calculations of the cross section [75-77]. The WW cross section is calculated 
at NLO accuracy in QCD using MCFM. The cross section for non-resonant gluon-fusion production is 
calculated at LO accuracy with GG2WW, including both WW and ZZ production and their interference. 

Top-quark event generation uses Powheg-i-Pythia6, except for the single-top t-channel process tqb, for 
which AcerMC [78]-i-Pythia6 is used. The WW background is also modelled using Powheg-i-Pythia6. For 
WW, WZ and ZZ backgrounds with two additional jets produced, the Sherpa generator is used for event 
modelling. The W{Zly*) process is simulated with Sherpa and Powheg-i-Pythia8, with tUy* extending 
down to the kinematic threshold and lepton masses included in the modeling of the y* decay. The Wy and 
Drell-Yan processes are modelled using Alpgen-i-Herwig with merged LO matrix element calculations 
of up to five jets. The merged samples are normalised to the NLO calculation of MCFM (for Wy) or 
the NNLO calculation of DYNNLO [79, 80] (for Zjy*). A Sherpa sample is used to model the Zy 
i£y background. The cross section of this process is normalised to NLO using MCFM. The VT-i-jets 
background shape and normalisation are derived from data, as described in Section 5.2, so no simulated 
VT-P jets events are used. 
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Table 2: Monte Carlo generators used to model the background processes in the H —» WW —» Mv analysis. All 
leptonic decay branching ratios (e, t) of the W and Z bosons are included in the product of cross section (cr) and 
branching ratio (BR). 


Background 

MC generator 

cr-BR(pb) 

tt 

Powheg -i- Pythia6 

26.6 

tW 

Powheg -i- Pythia6 

2.35 

tqb 

AcerMC -I- Pythia6 

28.4 

tb 

Powheg -i- Pythia6 

1.82 

qqlg WIT 

Powheg -i- Pythia6 

5.68 

gg WW 

GG2WW-P Herwig 

0.20 

QCD 1T1T-P2jets 

Sherpa 

0.568 

EW ITIT-P2jets 

Sherpa 

0.039 

Z/y*-f-jets (mu > 10 GeV) 

Alpgen -I- Herwig 

16.5 X 10^ 

EW Z/y* (includes f-channel) 

Sherpa 

5.36 

2(*)2(*) ^ 4t 

Powheg -i- Pythia8 

0.73 

WiZir)(m^z/r)<TGey) 

Sherpa 

12.2 

Zy{p\ > 1 GeV) 

Sherpa 

163 

Wy 

Alpgen -i- Herwig 

369 

Higgs boson {mu = 125 GeV) 

Powheg -i- Pythia8 

0.60 


4.2.2 Background processes for the H -» WW -»tvqq analysis 

Several different Monte Carlo generators are used to simulate the background to the H WW —> tvqq 
process. The processes used to model the background in the analysis are shown in Table 3. In general, the 
treatment follows that of the H WW —> tvtv analysis, with the exceptions described here. 

The VT+jets background is modelled with the Sherpa generator version 1.4.1. In order to have enough 
events for a background prediction at high mass, the Sherpa samples are generated in multiple bins of 
Pj. The bin boundaries are: 40-70 GeV, 70-140 GeV, 140-280 GeV, 280-500 GeV, and > 500 GeV. An 
inclusive sample is used for < 40 GeV. Samples of W bosons with only electroweak vertices are also 
generated to ensure sufficiently good modelling of this background in the VBF topology. 

The top-quark background is modelled using the same generators as in the H WW —> tvtv analysis. 
Events in the tt sample are reweighted according to the pj of the tt system and the individual top quarks 
to improve the kinematic agreement between the data and the Powheg prediction, following the prescrip¬ 
tion outlined in Ref. [81] based on the measurements of Ref. [82]. This treatment is not needed for the 
H WW —> tvtv analysis because the distributions affected are primarily the number of jets and the 
jet pj, and the analysis is not sensitive to either of these because of the normalisation of the top-quark 
background individually in each jet bin. The Z-i-jets background is also generated via Sherpa and, like the 
Sherpa 1T-p jets background, uses samples binned in p^, with a binning identical to the p^ used for the 
IT4-jets samples. 
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Table 3: Monte Carlo generators used to model the background processes in the H —» WW —> {vqq analysis, and 
the associated cross sections cr. Leptonic decay branching ratios BR of the W and Z bosons are not included in the 
number quoted unless explicitly indicated in the process name. 


Background 

MC generator 

cr • BR (pb) 

tt 

PowHEG -1- Pythia6 

252.9 

tW 

PowHEG -1- Pythia6 

22.4 

tqb 

AcerMC -i- Pythia6 

28.4 

tb 

PowHEG -1- Pythia6 

1.82 

W ^ €v 

Sherpa 

35.6 X 10^ 

W ^ iv VBF 

Sherpa 

12.6 

Wy —> bvy 

Alpgen -I- Herwig 

369 

Z ^ U 

Sherpa 

3.62 X 10^ 

Zy Uy{p\ > 10 GeV) 

Sherpa 

96.9 

IFW 

Herwig 

32.5 

WZ 

Herwig 

12.0 

ZZ 

Herwig 

4.69 


The Herwig generator is used for the WW, WZ, and ZZ processes. These samples are produced with 
inclusive vector boson decays and a single-lepton filter at the event generation stage. 


5 The H WW ^ tvtv analysis 

In the H WW —> channel, the final sfafe is fwo oppositely charged lepfons and fwo neufrinos, which 

are reconsfrucfed as missing fransverse momenfum. Additional jels may be presenl from QCD radiation or 
from Ihe scaflering quarks in fhe VBF produclion mode. The analysis described here is similar fo Ihe one 
designed lo sludy Ihe Higgs boson wilh mn ~ 125 GeV in Ihe WW^Mv final slate [83], wilh adaplafions 
made lo enhance Ihe sensifivily for a high-mass Higgs boson. 


5.1 Event selection 

The evenl is required lo have Iwo oppositely charged leptons and no additional lepton wilh pj > 10 GeV, 
wilh Ihe higher- and lower-px leptons respectively satisfying pj > 22 GeV and pj > 10 GeV. Bolh leptons 
musl satisfy Ihe qualify criteria discussed in Section 3. Background from low-mass resonances conslilules 
a significanl conlribulion, and is rejected by requiring > 10 GeV in Ihe same-flavour channel and 
mu > 12 GeV in Ihe differenl-fiavour channel, in which resonances decaying to tt may conlribule. In Ihe 
same-flavour channel, a veto on Z bosons is applied by requiring \mu - mz\ >15 GeV. These criteria form 
Ihe preselection. 
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Table 4: Event selection criteria used to define the signal regions in the H —> WW —» Mv analysis. The criteria 
specific to different-flavour (DF) and same-flavour (SF) channels are noted as such; otherwise, they apply to both. 
Preselection applies to all A^jet categories. In the >2 jets category, the rapidity gap is the rapidity range spanned by 
the two leading jets. 


Category 

Yet=0 

Ajet = 1 

IV 

Preselection 

Two isolated leptons {£ = e,p) with opposite charge 
> 22 GeV, > 10 GeV 

DF: m££ > 10 GeV 

SF: m££ > 12 GeV, | m££ - mz \ > 15 GeV 

Lepton pj 

p“>60GeV 

p™blead>30GeV 

p“>55 GeV 

p^ublead > 35 QgY 

p“ > 45 GeV 
p™bicad>20GeV 

Missing transverse 
momentum 

DF: p™"* > 45 GeV 
SF: > 45 GeV 

SF: > 65 GeV 

DF: p™"" > 35 GeV 
SF:£™“>45GeV 
SF: p™“ > 70 GeV 

DF:£™“ >25 GeV 

T.calo 

SF:£™;jl^>45GeV 

General selection 

p">60GeV 

' 

II 

o 

Nb-ia = 0 
p“<40GeV 


rtijj > 500 GeV 
Ayjj >4.0 

No jet (pj > 20 GeV) in rapidity gap 
Both £ in rapidity gap 


H WW Mv m{f > 60 GeV m££ > 65 GeV DF: m{£ > 60 GeV, SF: m££ > 45 GeV 

topology I^p££ < 1.35 < 1.35 ls.p££ < 1.85 


The signal and background compositions depend strongly on the final-state jet multiplicity (A^jet)- For 
A^jet = 0, the signal is predominantly from the ggF process, and WW events dominate the background. 
For A^jet = 1, both the ggF and VBF signal processes contribute, and the large majority of background 
events are from WW and top-quark events, which contribute approximately equally to the background. For 
A^jet > 2, the signal originates mostly from the VBF process and top-quark events dominate the background. 
The analysis is consequently divided into Vjet = 0, 1 and > 2 categories. 

The event selection in the various jet multiplicity categories is optimised using the BumpHunter [84] pro¬ 
gram, maximising the quantity s/ A/(FT(AAj^, where s and b are the numbers of signal and background 
events, respectively, and Ab represents the systematic uncertainty on the background. The value Ab = 10% 
is used. The optimisation has also been performed with Ab = 20% to test for sensitivity to the assumed 
systematic uncertainties, but the resulting selection is not significantly different from the one adopted. The 
optimisation is performed separately for the different- and same-flavour channels. The optimised event 
selection criteria that define the signal regions (SRs) in the analysis are summarised in Table 4. 

Owing to the topology of // —> WW Mv events, a selection on the missing transverse momentum 
is useful. In the different-flavour channel in both the Vjet = 0 and Vjet = 1 categories, requirements are 
imposed on In the same-flavour channel in these Vjet categories, selections on and A'™®® are 

used since, as explained in Section 3, these quantities efflciently reject Drell-Yan events. In the Vjet > 2 
category, A"™thresholds are used in both the different- and same-flavour channels. Selection using 
or in this category rejects a large fraction of signal events and is not optimal; they are therefore not 
used. 
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In the A^jet = 0 category, additional requirements on the pj of the dilepton system p^^and on mu are applied. 
In the A^jet = I category, a A-jet veto is applied to suppress the top background, and a selection on mu is 
imposed. To orthogonalise the A^jet = 0 and A^jet = 1 signal regions with respect to the WW control regions 
(Section 5.2), the pseudorapidity difference Arju between the two leptons is required to be smaller than 
1.35. 

The A^jet > 2 category is optimised to extract the Higgs boson signal produced via vector-boson fusion. The 
invariant mass mjj of the two highest-px jets, referred to as the tagging jets, is required to be larger than 
500 GeV. The magnitude of the rapidity difference between the tagging jets, Ayjj , is required to be larger 
than 4.0. In addition, the event must have no additional jets with px > 20 GeV within the rapidity gap 
of the tagging jets, while both leptons are required to be within this rapidity gap. A A-jet veto is applied, 
and the total transverse momentum p™ in the event is required to be smaller than 40 GeV. The quantity 
p“ is defined as the magnitude of Py1] p^^'^ where the sum is over all jets that pass the 
nominal analysis jet selection. Selections on mu are applied as in the Ajet < 1 categories, and Arju < 1-85 
is required. For a Higgs boson with mu - 300 GeV and the ratio of ggF and VBF cross sections predicted 
by the SM, 83% of the total signal selected in the Vjet > 2 category is produced by the VBF process. In the 
Vjet = 0 and Vjet = I categories, these fractions are 2% and 12%, respectively. The signal fractions from 
the VBF process increase with increasing mu- 

The discriminant used to derive the final results in this analysis is the transverse mass nix, defined as: 

nix = - Ipf + (2) 


where £■“ = yjlPjp- + ni^^. 

5.2 Background determination 

The major backgrounds in this analysis are top-quark and WW production, with additional contributions 
from W/Z+jets, multijets, and the diboson processes WZ, Wy, Wy*, and ZZ. The top-quark and WW 
backgrounds are normalised to data in control regions (CRs) defined by criteria similar to those used 
for the SR, but with some requirements loosened or reversed to obtain signal-depleted samples enriched 
in the relevant backgrounds. This normalisation is done through a simultaneous fit to the signal region 
and all control regions, as described in Section 9.1. This fit uses the complete background prediction in 
each region in order to account for the presence of other backgrounds and the potential small presence 
of signal. In particular, any background whose normalisation is determined by a control region is scaled 
by the same normalisation factor in all signal and control regions, not just its own control region. The 
following subsections describe the methods used to estimate the most important backgrounds, namely, 
WW, top-quark events, and IT-i-jets, in more detail. The Drell-Yan and non-ITIT diboson backgrounds 
are small, and their predictions are computed from simulation. The small background from the Higgs 
boson with mu ~ 125 GeV is also included. The predicted cross section, branching ratio, and kinematics 
for the SM Higgs boson are used. With few exceptions, the background estimates use the same techniques 
as Ref. [83]. They are described there in more detail, and summarized here. 
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Table 5: Event selection criteria for the Njf-i — 0 and — 1 WW control regions in the H —» WW —» Mv analysis. 
The criteria that are different with respect to the SR definition are shown. Only the different-flavour final state is 
used. 


Category 

Ajet^O 

Ajet = 1 

Lepton transverse momentum 

22 GeV 

^sutlead > ^5 QgY 

Missing transverse momentum 

> 20 GeV 

> 35 GeV 

General selection and 

pf >35GeV 

iVi-jet=0 

H WW Mv topology 

mu >75 GeV 

mu >75 GeV 


Arju > 1.35 

Ariu> 1.35 


5.2.1 WW background 

In the A/jet < 1 categories, the WW background is normalised using a CR defined with the selection sum¬ 
marised in Table 5. To orthogonalise the WW CRs to the A^jet = 0 and A^jet = 1 SRs, the selection on Aijfe 
is reversed with respect to the SR definitions: Arjge > 1.35 is required. Only the different-flavour final 
states are used to determine the WW background, and the purity is 70.5% and 40.6% in the TTjet - 0 and 
A/jet = 1 categories, respectively. The normalisation factors obtained from the simultaneous fit to the signal 
and control regions are 1.18 ± 0.04 for the Njst -0 CR and 1.13 ± 0.08 for the Njst - 1 CR, where the 
uncertainty quoted includes only the statistical contribution. The high normalisation factor for WW events 
with zero jets has been studied in Ref. [83], and results from poor modelling of the jet veto efflciency. 
The WW prediction in the A^jet > 2 category is taken from simulation, because it is difficult to isolate a 
kinematic region with a sufficient number of WW events and a small contamination from the top-quark 
background. 

Figure 1 shows the mj distributions in the A^jet < 1 WW CRs. Normalisation factors obtained from the top 
CRs as well as from the WW CRs have been applied to these distributions. 


5.2.2 tt and single top background 

Top-quark events can be produced as a tt pair, or in association with a W boson or another flavour of quark. 
In the H WW —> Mv analysis, contributions from tt and single-top events are estimated together, 
with their relative contributions determined by the predicted cross sections and MC simulation of the 
acceptances, since it is not easy to kinematically separate the two processes and the contribution from 
single top is relatively small. 

Owing to the difficulty of defining reasonably pure control regions in the Njet - 0 category, the top-quark 
background in this category is not estimated from the likelihood fit. The jet veto survival probability 
(JVSP) procedure, described in more detail in Ref. [85], is employed instead. In this method, the normal¬ 
isation is derived from the top-quark event yield determined in a control region defined by events with a 
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Figure 1: Transverse mass distribution in the A^jet = 0 (left) and A^jet = 1 (right) WW control regions of the 
H —» WW —» Mv analysis. Only the different-flavour final state is used. In each hgure, the last bin contains 
the overflow. The combined statistical and systematic uncertainties on the prediction are shown by the hatched band 
in the upper pane and the shaded band in the lower pane. Normalisation factors obtained from a comparison of data 
and prediction have been applied in these figures. 


different-flavour opposite-sign lepton pair, any number of jets, and p™*® >45 GeV. This sample is domin¬ 
ated by top-quark events. The estimated top-quark event yield is the total number of events VcR passing 
this selection minus the expected contribution Bqr from other processes. The theoretical cross sections 
and acceptances from MC simulation are used to calculate Bqr, except the IV-i-jets background, for which 
the data-derived estimate described later in this section is used. The resulting estimated top-quark event 
yield is multiplied by the fraction cq of top-quark events with no reconstructed jets obtained from simu¬ 
lation in the CR. This fraction is corrected using data from a second CR defined like the first, with the 
additional requirement of at least one ti-tagged jet. The fraction of events in this CR with zero jets in 
addition to the ti-tagged one is measured in both data and simulated top-quark events, denoted and 
respectively. Using these inputs, the estimated number of top-quark background events in the 
Vjet = 0 signal region is estimated as: 

tVtop = (A^cr - Bcr) • eo • • Crest, (3) 

where Crest is the efficiency of the Vjet - 0 selection requirements applied after the jet veto, derived from 
simulated top-quark events. The theoretical uncertainties on the quantities derived from top-quark MC 
simulation, namely cq, and Crest, are described in Section 7. 

In the Vjet = 1 and Vjet > 2 categories, the normalisation of the top-quark background is determined from 
control regions. As with the WW CR, and unlike the Vjet - 0 CRs, these are included in the simultaneous 
fit with the signal regions. These CRs are defined identically to the respective signal regions, except 
that the threshold is lowered to 20 GeV and the veto on ti-tagged jets is inverted to require exactly 
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Figure 2: Transverse mass distribution in the A^jet = 1 (left) and A^jet > 2 (right) top control regions of the 
H —> WW Mv analysis. In the A^jet = 1 category only the different-flavour final state is used; in the A^jet > 2 
category different-flavour and same-flavour final states are used. In each figure, the last bin contains the overflow. 
The combined statistical and systematic uncertainties on the prediction are shown by the hatched band in the up¬ 
per pane and the shaded band in the lower pane. Normalisation factors obtained from a comparison of data and 
prediction have been applied in these figures. 


one Zt-tagged jet with pj > 25 GeV. The purity is 96.5% in the A^jet - 1 category and 90.7% in the 
TTjet > 2 category. In the Vjet = 1 category, only the different-flavour final states are used to obtain the 
normalisation. In the Vjet > 2 category same-flavour and different-flavour final states are used to increase 
the number of events and thereby improve the statistical precision. The normalisation factors obtained 
from the simultaneous fit to the signal and control regions are 1.05+0.03 for the Ajet - 1 CR and 0.92+0.06 
for the Vjet > 2 CR, where the uncertainty quoted includes only the statistical contribution. Figure 2 shows 
the niT distributions in the Vjet = 1 and Vjet > 2 top CRs. The normalisation factors have been applied in 
these distributions. 


5.2.3 W+jets and multijet background 

The procedures to estimate the VF+jets and multijet backgrounds using data are described in more detail in 
Ref. [83] and summarised here. The 1V+jets background contribution is estimated using a control sample 
of events in which one of the two lepton candidates satisfies the identification and isolation criteria used 
to define the signal sample (these lepton candidates are denoted “fully identified”), and the other lepton 
fails to meet these criteria and satisfies a less restrictive selection (denoted “anti-identified”). Events in 
this sample are otherwise required to satisfy all of the signal selection criteria. The dominant component 
of this sample (85% to 90%) is 1V+jets events in which hadronic activity produces an object reconstructed 
as an anti-identified lepton. It may be either a non-prompt lepton from the decay of a hadron containing a 
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heavy quark, or a particle from a jet reconstructed as a lepton candidate. 

The VT+jets contamination in the SR is determined by scaling the number of events in the control sample 
by an extrapolation factor, which is measured in a data sample of jets produced in association with Z 
bosons reconstructed in either the e^e~ or the final state. Kinematic vetoes reduce contamination 
from ZZ and ITZ events, and the expected remaining contribution is subtracted. The extrapolation factor 
is the ratio of the number of fully identified leptons to fhe number of anfi-idenfified lepfons, measured 
in bins of anfi-idenfified lepfon pj and rj. To accounf for differences befween fhe jefs associafed wifh W 
and Z boson producfion, fhe exfrapolafion facfors are measured in simulafed IT-i-jefs and Z-t jefs evenfs, 
and fhe ratio of fhe fwo exfrapolafion facfors is mulfiplied by fhe one measured in fhe Z-i-jefs dafa as a 
correction. The cenfral value of fhe correcfion facfor is close fo unify; differences among Monfe Carlo 
generafors for fhis ratio of abouf 20% are observed and are faken as a sysfemafic uncerfainfy. 

The background in fhe SR due fo mulfijefs is defermined using a confrol sample fhaf has fwo anfi-idenfified 
lepfon candidafes, buf ofherwise safisfies all of fhe SR selection criferia. An exfrapolafion facfor is esfim- 
afed using a mulfijef sample in dafa and applied fwice fo fhe confrol sample. The sample used fo defermine 
fhe exfrapolafion facfor is expecfed fo have a similar sample composifion in ferms of heavy-flavour jefs, 
lighf-quark jefs and gluon jefs as fhe confrol sample. Since fhe presence of one misidenfified lepfon in a 
mulfijef sample can change fhe sample composifion wifh respecf fo a mulfijef sample wifh no lepfon se- 
lecfion imposed, correcfions fo fhe exfrapolafion facfor are made fhaf fake info accounf such correlafions. 
These are evaluafed using a mulfijef Monfe Carlo sample and range from 1.0 to 4.5 depending on fhe 
lepfon pt and flavour. The uncerfainfy on fhese is 30-50%, wifh fhe dominanf confribufion being from fhe 
heavy-quark cross sections. 


6 The H WW tvqq analysis 

In fhe H WW —> ivqq channel fhe final sfafe consisfs of one W boson decaying info a quark-anfiquark 
pair leading fo a pair of jefs, wifh fhe ofher W boson decaying info a charged lepfon and a neufrino 
(W (v, wifh i - e or p). This channel is parficularly sensifive in searching for a Higgs boson wifh a 
mass greafer fhan fwice fhe W boson mass since tuh can be reconsfrucfed on an evenf-by-evenf basis and 
used as fhe discriminanf fo search for a signal. This evenf-by-evenf reconsfrucfion is done using kinemafic 
consfrainfs fhaf provide an esfimafe of fhe componenf of fhe neufrino momenfum along fhe beam axis and 
require signal jefs in fhe evenf fo be consisfenf wifh coming from a hadronic W decay. 


6.1 Event preselection and categorisation 

Evenfs are required fo have exacfly one reconsfrucfed lepfon candidafe {e or p) wifh pT > 25 GeV; no 
additional lepfon wifh pj > 15 GeV is allowed. The selecfed lepfon musf mafch fhe objecf fhaf friggered 
fhe evenf. Evenfs in fhe SR are required fo have > 60GeV in order fo suppress mulfijef processes 

while refaining a high signal efficiency. 

Jefs are used fo distinguish befween ggE and VBE producfion as well as fo reconsfrucf fhe hadronic W 
boson decay. Anfi-k( jefs are selecfed wifh pj > 30 GeV and \q\ < 4.5, and large-/? jefs are selecfed wifh 
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Pt: > 100 GeV, \r]\ < 1.2, and mj > 40GeV, where the J subseript indieates a large-/? jet. Both the anti-^, 
and large-/? jets are required to be separated from the eharged lepton by A/? > 0.3. There is no explieit 
overlap removal between anti-^f and large-/? jets. 

The momentum of W bosons from the Higgs boson deeay inereases with inereasing Higgs boson mass. 
This feature leads to a progressively smaller opening angle between the jets produeed by the W boson 
deeay, making the jets diffieult to distinguish using standard jet reeonstruetion algorithms. To mitigate the 
resulting loss in signal effieieney, the hadronie W deeay may be reeonstrueted from either two anti-k? jets 
or one large-/? jet eonsistent with originating from a W boson deeay. 

In hadronie W boson deeays reeonstrueted from two anti-k^ jets, the best eandidate jet pair is referred to 
as the “W jets”. The two jets with an invariant mass elosest to the W boson mass are taken to be the W 
jets, unless there is more than one jet pair with |m,y - mw| < 15 GeV. In that ease, the pair having the 
highest Pt is ehosen. Categorisation by produetion mode is done prior to identifieation of the hadronie W 
boson, so jets identified as the VBF tagging jets aeeording to the proeedure deseribed below are exeluded. 
Hadronie W boson deeays are identified using a single large-/? jef if fhere is one wifh pj > 100 GeV and a 
mass eloser fo fhe W boson mass fhan fhe invarianf mass of fhe besf dijef pair. In fhis ease, fhe large-/? jef 
replaees fhe W jefs as fhe eandidafe for fhe hadronieally deeaying W boson. 

Evenfs are elassified info fwo eafegories designed fo disfinguish befween fhe ggF and VBF produefion 
modes, based on fhe number of jefs in fhe evenf and fhe properfies of fhose jefs. In fhe eafegory designed 
fo be sensitive fo fhe ggF produefion mode, referred fo as fhe ggF seleefion, all evenfs are required fo have 
af leasf fwo anfi-k? jefs or af leasf one large-/? jef, and fail fhe VBF seleefion. 

In fhe seeond eafegory, designed fo be sensifive fo fhe VBF produefion mode, evenfs are required fo have 
af leasf four anfi-k( jefs or af leasf fwo anfi-k, jefs and one large-/? jef. Orfhogonalify befween fhe ggF and 
VBF eafegories is ensured by idenfifying fhe fwo anfi-k; jefs ji and 72 wifh fhe largesf invarianf mass, and 
assigning fhe evenf fo fhe VBF (ggF) eafegory if fhese jefs pass (fail) fo meef eriferia, referred fo as fhe 
VBF seleefion, eharaeferisfie of fhe forward jefs produeed by fhe VBF proeess. This VBF fagging jef pair 
is required fo have an invarianf mass mj^J^ > 600 GeV, wifh fhe leading jef pi > 40 GeV, and be well 
separated in rapidity sueh fhaf - If/;, - > 3. If fhe ratio of fhe ggF and VBF eross seefions 

is as prediefed by fhe SM, 63% of signal evenfs passing fhe full VBF preseleefion are produeed via VBF, 
and 93% of signal evenfs passing fhe ggF preseleefion are produeed via ggF. 

Vetoes, based on fhe presenee of k-jefs in fhe evenf, rejeef tt baekground. If bofh of fhe Wjefs are k-fagged, 
fhe evenf is vefoed. If only one of fhe W jefs is k-fagged, fhe even! is kepf to mainfain signal effieieney 
sinee a large fraelion of jefs from W ^ cs deeays are k-fagged. If any ofher jef in fhe evenf is k-fagged, 
ineluding fhe VBF lagging jefs, fhe even! is vefoed. If a large-/? jef is used to reeonsfruel fhe W boson, 
evenfs wifh k-fagged jefs oulside of A/? = 0.4 from fhe axis of fhe large-/? jef are vefoed. No flavour 
lagging is applied to large-/? jefs. 

Furlher seleelions are applied fo ggF and VBF seleeled evenfs. In bofh eafegories, eaeh of fhe W jefs is 
required fo have \r]\ < 2.4 and Iheir invarianf mass fo be in fhe range 65 GeV < nijj/j < 96 GeV, fhaf 
is, elose fo fhe W boson mass. Additionally, for hadronie W boson eandidales reeonslrueled from fwo 
anli-k, jefs, one of fhe fwo W jefs is required fo have pj > 60 GeV in bofh ggF and VBF seleeled evenfs. 
Furlher requiremenls are imposed on fhe azimulhal separalion of reeonslrueled objeels whieh exploil fhe 
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Table 6: Summary of event preselection in the H —> WW —> £vqq analysis. The “tagging jets” yj and j 2 are the pair 
of anti-A:, jets with the highest invariant mass among all pairs in the event, and yj is the higher-px jet in the pair. The 
decay topology selection differs between events in which the hadronic W boson candidate is reconstructed as a pair 


of jets jj or as a single large-/? jet J. For the jet-pair topology, the leading jet is denoted yjead and if only a single jet 
j is referenced, the requirement is applied to both jets. 

Object selection 

1 isolated charged lepton {e or p): px > 25 GeV, I 77 I < 2.4 

> 60GeV 

jet: px > 30 GeV, \t]\ < 4.5 
large-/? jet: px > 100 GeV, Ipl <1.2 

VBF selection 

(> 4 jets) or (> 2 jets - 1 - > 1 large-/? jets) 

> 600 GeV 
p{‘ > 40 GeV 

ApO'i,/ 2 ) > 3.0 

ggF selection 

not VBF tagged and (> 2 jets or > 1 large-/? jet) 

Further selection, hadronic 

W boson reconstructed as: jet pair large-/? jet 

Decay topology 

p^^** > 60GeV 

A0(77) <2.5 

A(f>{jJ)> \.0 A(/){J,£)> l.O 

A 0 aE™-o)>l-O A0(7,E™;V>1-O 

A0(AESo)<2.5 

/^-tagging 

veto events with: 

both W candidate jets /^-tagged /^-tagged jet with AR{j, J) > 0.4 
or any other jet / 7 -tagged 

VF-mass window 

65 GeV < Mjj < 96 GeV 65 GeV <mj< 96 GeV 


decay topology of signal events to improve the expected sensitivity. A summary of the event preselection 
is shown in Table 6. 

The signal region is subdivided into exclusive categories which separate sources of signal and background. 
In addition to the ggF and VBF selection which separates the two signal production modes, the signal 
regions are separated by the flavour of the charged lepton and the sign of its charge. Electrons and muons 
are affected differently by the multijet background, and the categories with positively-charged leptons 
have a higher proportion of IF-i-jets background because of the charge asymmetry of W production in pp 
collisions. 
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6.2 WW invariant mass reconstruction 


The invariant mass of the WW system is reconstructed from the four-momenta of the two W boson can¬ 
didates. The reconstructed invariant mass is denoted niiyjj regardless of whether the hadronic W boson is 
reconstructed from a jet pair or a large-7? jet. The reconstruction of the leptonic W boson decay relies on the 
charged lepton and neutrino four-momenta. The complete four-vector is measured for the charged lepton, 
and the provides the transverse components of the neutrino momentum. The neutrino longitudinal 
momentum is computed using the quadratic equation resulting from the mass constraint m{{v) = m{W). 
In the case of two real solutions of this equation, the solution with the smaller |p^| is taken. In the case 
of complex solutions, only the real part of the solution is taken. Based on signal simulation, this proced¬ 
ure has been shown to give the correct p^ solution in 60-70% of events after the preselection, depending 
on the event category and the Higgs boson mass. The experimental mass resolution of the reconstructed 
Higgs boson varies from ~ 30 GeV for mu = 300 GeV to ~ 60 GeV for mn - 1 TeV. For mn - 420 GeV, 
the mass resolution is about the same as the width of a SM Higgs boson at that mass, ~ 36 GeV. 


6.3 Signal region selection 

The sensitivity to a heavy Higgs boson in the H —> WW —> ivqq channel is improved by applying event 
selection in addition to the preselection described in Section 6.1, as a function of the Higgs boson mass 
hypothesis. 

The mass-dependent optimised selection is based on a set of kinematic quantities that discriminate between 
the signal from a hypothetical CP-even scalar and the background. For the ggF selection, these are the 
leading jet pj, the subleading jet px, the large-/? jet px, the lepton px, ^<piv and For the VBF 
selection, these are Acpjj, and px balance, where px balance is defined as 


PT.balance — 


(p^ +P'' + + pJ'*)^ 

+ p^ + p-!^ + p!^ 


(4) 


with 73 and 74 representing the VT-jets, to distinguish them from the tagging jets 71 and 72 . In the case of 
large-/? jet events, the terms representing 73 and 74 are replaced with a single term p;[ which represents the 
large-/? jet momentum instead. Also, the Acl)jj selection is only applied to events in which two resolved 
jets form the hadronically decaying W boson. Fewer criteria are used for the VBF selection than the ggF 
selection because of the smaller event yields in the VBF channel. 

The selection is optimised as a function of the Higgs boson mass through a two-step procedure. In the 
first step, the selection that optimises the expected signal significance are found in 100 GeV incremenfs of 
fhe mass hypofhesis mn- The expected signal significance is defined as s/ + b, where 5 is fhe number 
of expecfed signal evenfs and b is fhe number of expected background evenfs. Ofher esfimafors for fhe 
significance have been tested {s/ and s/ '\lb + Ab wifh Ab = 10% and 30%) and shown fo provide 
fhe same opfimal selection. Since fhe majorify of signal evenfs are localised fo a region in niiyjj thaf is 
small compared wifh fhe overall fif region, fhe significance calculafion does nof include all signal and 
background evenfs, buf rafher only evenfs in which niiyjj is wifhin a specified range around mh, defined 
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as the region that contains 90% of the signal events. The resulting selection criteria become stricter 
with increasing mn as the decay products are produced at higher momenta, allowing greater background 
rejection while maintaining good signal efficiency. 

In the second step, the selection is slightly relaxed, because the optimal selection, particularly for low mu, 
typically causes the peak of the signal mcvjj distribution to coincide with that of the expected background, 
thus reducing the sensitivity of the analysis because of large systematic uncertainties in describing the 
turn-over region in mivjj. Typically, the optimal value for the next lower mu increment is used. 

6.4 Background estimation 

In all signal regions, the background is expected to be dominated by Vk-i-jets production, with other import¬ 
ant contributions from tt, single top, and multijet production that can be selected owing to the presence 
of leptons from heavy-flavour decays or jets misidentified as leptons. Diboson events, including WW, 
VTZ, ZZ, Wy, and Zy, as well as Z-i-jets events, contribute at a smaller level and are also accounted for as 
backgrounds. 

The Vk-i-jets and top-quark backgrounds are modelled using simulation but their normalisations (one for 
VT-rjets and one for top-quark backgrounds) are determined through a simultaneous fit to the signal and 
control regions, similarly to what is done for the £v{v final state. The profile likelihood fit is described in 
more detail in Section 9.1. Multijet backgrounds are estimated using a CR selected in data to be enriched 
in leptons produced by hadronic activity. The small additional background from Z-i-jets and dibosons, 
including Zy and Wy, are estimated using their theoretical cross sections with simulation for the event 
selection acceptance and efficiency. 

6.4.1 W+jets and top-quark background 

IT-i-jets and top-quark production are the most important backgrounds in the H WW —> ivqq analysis. 
Their normalisations are set, and the mcyjj shape corrected, using the data observed in the corresponding 
CRs, defined below. 

The W control region (WCR) is defined similarly to the SR, but with the signal contributions suppressed 
by rejecting events with a dijet mass, or a large-R jet mass, consistent with the hadronic decay of a Vk 
boson. The WCR for the ggF selection is defined by the upper and lower sidebands to the reconstructed 
W boson mass. 


52 GeV < mjj < 65 GeV, ggF lower sideband 
96GeV < mjj < 126 GeV, ggF upper sideband 


(5) 

( 6 ) 


using the dijet or large-/? jet mass closest to the W mass, as described in Section 6.1. The corresponding 
sidebands in the WCR for the VBF selection are 


43 GeV < mjj < 65 GeV, VBF lower sideband 
96 GeV < mjj < 200 GeV, VBF upper sideband. 


(V) 

( 8 ) 
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The width of the sidebands is increased in the VBF selection compared to the ggF selection to improve the 
statistical precision of the background estimate. For both the ggF and VBF selection, the number of back¬ 
ground events is similar in the upper and lower sidebands, so no statistical bias is generated by using them 
jointly to define a common background normalisation. Separate CRs are defined for the ggF and VBF 
selection, but all lepton flavours and charges as well as the large-/? and dijet W reconstruction topologies 
are merged into a single CR. The CR selection also follows the mass-dependent SR selection described in 
Section 6.3. The purity of the WCR depends on the particular selection, but varies between about 70% and 
over 80%, with higher purities for higher mh selection and for ggF compared to VBF. The VT-i-jets norm¬ 
alisation factors are consistent with unity and stable with respect to the m//-dependent selection and ggF 
vs. VBF selection. The normalisation factors resulting from the simultaneous fit to the signal and control 
regions range from 0.9 ± 0.1 to 1.3 ± 0.3. The quoted uncertainties include statistical and experimental 
systematic components. The highest values are found with the VBF and high-mn selection. 

Figure 3 shows the m^yjj distributions observed in the ggF and VBF WCRs. Non-lT-boson contributions 
are subtracted from the data, and the resulting distribution is compared to the prediction from simulated 
VT-rjets events, after normalising the prediction to the data. The observed distributions differ substan¬ 
tially from those predicted in both the ggF and VBF WCRs. An alternative MC sample, generated using 
Alpgen -I- PythiaO, results in better agreement with the data in the WCR, but does not have enough events 
for a statistically precise prediction in the signal region for large values of nin- 

To correct the observed mismodelling, simulated W-i-jets events are reweighted using the rtifyjj distribu¬ 
tion observed in the WCR. In order to obtain a smooth function to use for the reweighting, the ratio of 
the data to the prediction is fit with polynomial functions. The degree of the polynomial is chosen to have 
enough flexibility to yield a good fit quality, the fit range being restricted to m^yjj values where a statist¬ 
ically meaningful fit can be made. For the ggF WCR, shown on the left side of Figure 3, a second-order 
polynomial function is used and the fit is extended to meyjj = 1.7 TeV. Above that value of mfyjj, a constant 
function at the value of the polynomial at m^yjj =1.7 TeV is used. For the VBF WCR, shown on the right 
side of Figure 3, a third-order polynomial function is used, up to niiyjj = 0.9 TeV. Fits extending to higher 
values of nifyjj have been attempted but require either a more complex fitting function or have a visibly 
poor-quality fit to the data. For simplicity, a constant function is used for miyjj > 0.9 TeV, as illustrated in 
the figure. The value of the function is the value of the third-order polynomial at m^yjj = 0.9 TeV. 

The top-quark control region (TopCR) is designed to be as pure as achievable for the second largest 
background, tt WbWb €vjj + bb. The event topology of this background is similar to that of the 
Higgs boson signal, but contains two characteristic fj-jets. The TopCR is defined to be identical to the 
SR, but with the f^-jet veto reversed. As with the WCR, the TopCR region selection follows the mass- 
dependent signal region selection described in Section 6.3, and lepton flavours and hadronic W topologies 
are merged but separate TopCRs are defined for the ggF and VBF topologies. The purity of the TopCR 
is about 80% and does not depend strongly on the region-specific selection. Similarly, the value of the 
resulting normalisation factor is stable with respect to the kinematic selection and is consistent with unity 
within the uncertainties. The values of the normalisation factor found by the simultaneous fit to the signal 
and control regions range from 0.9 ± 0.1 to 1.3 ± 0.2. Both extremes occur in the VBF control region for 
m//-dependent selection for mu > 700 GeV, which is most subject to statistical fluctuations. 

The mtyjj distributions in the ggF and VBF TopCRs are shown in Figure 4. As in Figure 3, processes 
other than top-quark single and pair production are subtracted from the data, and the resulting distribution 
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Figure 3: Data and Monte Carlo comparison of the shape of the invariant mass of the WW system nifyjj in the 
ggF (left) and VBF (right) WCRs after the mn — 300 GeV selection for the H —» WW —» tvqq analysis. All the 
lepton flavour and charge categories are summed together. To isolate the effects of fT+jets background modelling, 
other contributions (top, diboson, Z+jets, multijet) are subtracted from the data. The Monte Carlo distributions are 
normalised to the remaining data. The ratio of the data to the Monte Carlo distribution is shown in the bottom panel, 
along with a red line showing the resulting weights that are applied to correct the Monte Carlo predictions in the rest 
of the analysis. 


is compared to the top-quark prediction after normalising the prediction to the data. Similarly to the WCR, 
differences in shape are observed, and simulated top-quark events are reweighted accordingly as a function 
of mcvjj. A first-order polynomial function is fit to the data. Since the purity of events with a top quark in 
the TopCR is very high, the fT-i-jets contribution in Figure 4 has no miyjj reweighting applied. 

For both the VF-i-jets and top-quark backgrounds, the event weights used to correct the MC simulated 
events are derived from corresponding CRs after the mn - 300 GeV selection, because this is the most in¬ 
clusive selection. The higher-mass criteria select subsets of the events accepted by the 300 GeV selection. 
The agreement between the data and MC distributions in the W and top-quark CRs when other mass- 
dependent selection criteria are applied is consistent with the results from the 300 GeV selection. The 
reweighting is event-by-event, using the fitted functions shown in Figures 3 and 4 and the m(yjj value of 
each simulated VF-i-jets or top-quark background event. The reweighting is applied after the preselection, 
and therefore propagates to all signal and control regions. For events with rriiyjj above the fitted range, the 
weight at the upper boundary of the fit is used. 

Half of the difference between the nominal and the reweighted miyjj distribution is taken as a systematic 
uncertainty on the m(yjj shape of these backgrounds. These uncertainties are included as Gaussian con¬ 
straints in the profile likelihood fit, which allows the fit to adjust the shape of these backgrounds using the 
shape of the data in the signal regions. 
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Figure 4: Data and Monte Carlo comparison of the shape of the invariant mass of the WW system nifyjj in the ggF 
(left) and VBF (right) TopCRs after the mn - 300 GeV selection for the H —» WW —> £vqq analysis. All the lepton 
flavour and charge categories are summed together. To isolate the effects of top-quark background modelling, other 
contributions (fk-njets, diboson, Z-i-jets, multijet) are subtracted from the data. The Monte Carlo distributions are 
normalised to the remaining data. The ratio of the data to the Monte Carlo distribution is shown in the bottom panel, 
along with a red line showing the resulting weights that are applied to correct the Monte Carlo predictions in the rest 
of the analysis. 


6.4.2 Multijet background 

The shapes of multijet background distributions are modelled using histograms derived from data samples 
selected similarly as for signal events, except that the lepton identification requirements are loosened, 
while the isolation requirement is not changed. In the electron channels, a loosened identification selection 
is applied to the data, with a veto for electrons selected using the standard criteria. In the muon channels, 
the impact parameter significance requiremenf is reversed. 

The normalisafion of the multijet background in a given event category is derived from a standalone 
template fit (separate from the final simulfaneous fif) fo fhe disfribufion wifhouf fhe requiremenf 
applied. The femplafe for fhe mulfijef background in each region is faken from dafa selecfed wifh modified 
lepfon selecfion as described above, buf ofherwise following the event selection of that signal or control 
region. The relative contributions of backgrounds other than the multijet background are fixed fo fheir SM 
expecfafions. 

The mulfijef background is relafively small, consfifufing befween 1% and 4% of fhe fofal expected back¬ 
ground depending on the selection used, with smaller contributions for the mass-dependent SR selection 
applied at higher values of nin- This background is also concentrated at low values of nifyjj, so its effect 
is reduced by the miyjj shape fit. 


26 



























































7 Systematic uncertainties 


This section describes the systematic uncertainties affecting the analysis results. Experimental and theor¬ 
etical uncertainties common to the H WW —> £v{v and H WW —> {vqq analyses are described first, 
followed by a discussion of uncertainties particular to each channel. 


7.1 Common experimental uncertainties 

The dominant sources of experimental uncertainties on the signal and background yields are the jet energy 
scale and resolution, and the f7-tagging efficiency. Ofher sources of uncertainfy include fhose on fhe scale 
and resolution of fhe lepfon energy or momenfum, lepfon idenfificafion and frigger efficiencies, fhe scale 
and resolufion of fhe missing fransverse momenfum, and fhe luminosify calibrafion. All experimenfal 
uncerfainfies are freafed by varying fhe objecf subjecf fo a parficular uncerfainfy and fhen re-running fhe 
full analysis. 

The Jef energy scale is defermined from a combinafion of fesf-beam dafa, simulafion, and in situ meas- 
uremenfs [27]. Ifs uncerfainfy is splif info several independenf cafegories: modelling and sfafisfical un- 
cerfainties on fhe exfrapolafion of fhe jef calibrafion from fhe cenfral region (q infercalibrafion), high-pj 
jef behaviour, Monfe Carlo non-closure uncerfainfies, uncerfainfies on fhe jef quark and gluon composi- 
fions and fheir calibrafions, fhe b-]&t energy scale uncerfainfies, uncerfainfies due fo modelling of in-fime 
and oul-of-fime pile-up, and uncerfainfies on in situ jef energy corrections. Some of fhese cafegories are 
furfher subdivided by fhe physical source of fhe uncerfainfy. For fhe anfi-k; jefs used in fhese analyses, 
fhe jef energy scale uncerfainfy ranges from 1% fo 7% depending on px and q. The resolufion varies 
from 5% fo 20%, and fhe relafive uncerfainfy on fhe resolufion ranges from 2% fo 40%. The lowesf- 
Pt jefs, immediafely above fhe jef selecfion fhresholds, have bofh fhe pooresf resolufion and fhe largesf 
uncerfainfy. 

The evaluafion of fhe b-jet tagging efficiency uses a sample dominafed by dileptonic decays of fop-quark 
pairs [34]. To improve fhe precision, fhis mefhod is combined wifh a second calibrafion mefhod based 
on samples confaining muons reconsfrucfed in fhe vicinify of fhe jef. The uncerfainties relafed fo ft-jef 
idenfificafion are decomposed info six uncorrelafed componenfs using an eigenvecfor mefhod [33], fhe 
number of componenfs being equal fo fhe number of px bins used in fhe calibrafion. The uncerfainfies 
range from < 1% fo 7.8%. The uncertainties on the misidentification rate for light-quark jets depend on 
Px and q, with a range of 9%-19%. The uncertainties on c-jets reconstructed as f7-jets range between 
6%-14% depending on the jet px. 

The reconstruction, identification, isolation, and trigger efficiencies for electrons and muons, as well as 
their momentum scales and resolutions, are estimated using Z—>ee,pp, >ee,pp, and lT^ev,py 
decays [22, 23, 25]. The uncertainties on the lepton identification and trigger efficiencies are smaller fhan 
1% excepf for fhe uncerfainfy on fhe elecfron idenfificafion efficiency, which varies befween 0.2% and 
2.7% depending on px and q. 

The changes in jef energy and lepfon momenfa due fo sysfemafic variafions are propagafed fo E™*®, such 
fhaf changes in fhe high-px objecf momenfa and in are fully correlafed. Addifional confribufions fo 
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the uncertainty arise from the modelling of low-energy particle measurements (“soft terms”) [35]. 
Calorimeter measurements of these particles, used in use calibrated clusters of cells not associated 

with reconstructed physics objects, with a noise threshold applied. The longitudinal and perpendicular 
components of the soft terms, defined with respect to the computed using hard objects, are smeared 
and rescaled to evaluate the associated uncertainties. The uncertainties are parameterised as a function of 
the magnitude of the vector sum pj of the high-px objects, and are evaluated in bins of the average number 
of interactions per bunch crossing. Differences of the mean and width of the soft term components between 
data and simulation result in variations on the mean of the longitudinal component of about 0.2 GeV. The 
uncertainty on the resolution of the longitudinal and perpendicular components is 2% on average. The 
systematic uncertainties related to the track-based soft term for are calculated by comparing the 
properties of p™'*® in Z —> ee,pp events in data and simulation as a function of the magnitude of the 
summed px of the leptons and jets in the event. The variations on the mean of the longitudinal component 
are in the range 0.3-1.4 GeV and the uncertainties on the resolution on the longitudinal and perpendicular 
components are in the range 1.5-3.3 GeV, where the lower and upper bounds correspond to the range of 
the sum of the hard px objects below 5 GeV and above 50 GeV, respectively. 

The uncertainty on the integrated luminosity is 2.8%. It is derived following the same methodology as in 
Ref. [86] from a calibration of the luminosity scale derived from beam-separation scans. 


7.2 Common theoretical uncertainties 

Theoretical uncertainties on the signal production cross section affect the H WW —»Mv and H —> WW —» {vqq 
analyses in the same way. These include uncertainties due to the choice of QCD renormalisation and fac¬ 
torisation scales, the PDF model used to evaluate the cross section and acceptance, and the underlying 
event and parton shower models. These are described and evaluated as a function of niH in Refs. [87, 88]. 

The QCD scale uncertainty on the inclusive signal cross sections is evaluated to be 8% for ggF and 1% 
for VBF production. The PDF uncertainty on the inclusive cross sections is 8% for ggF and 4% for VBF 
production. Uncertainties on the interference weighting of the CPS signal samples, described in Section 4, 
are also included. 


7.3 Uncertainties specific to the H WW (v(v analysis 

The uncertainties specific to the H WW —> Mv analysis arise primarily from the theoretical modelling 
of the signal acceptance in jet bins, the theoretical models used in the background predictions, and the 
additional consideration of uncertainties on the mx shape used in the likelihood fit. Statistical uncertainties 
on yields in control regions, and the effect of subtracting other processes from the control region yield, are 
also included in the total uncertainty on background yields predicted using control regions. 

Since the analysis is binned by jet multiplicity, large uncertainties from variations of QCD renormalisation 
and factorisation scales affect the predicted contribution of the ggF signal in the exclusive jet bins, and can 
cause event migration among bins. These uncertainties are estimated using the HNNLO program [79, 89] 
and the method reported in Ref. [90] for Higgs boson masses up to 1500 GeV. The sum in quadrature of 
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the inclusive jet bin uncertainties amounts to 38% in the 0-jet category and 42% in the 1-jet category for 
mu - 600 GeV. For mu - 1 TeV, these uncertainties are 55% and 46%, respectively. 

For the backgrounds normalised using control regions, theoretical uncertainties arise from the use of 
simulations of the background used in the extrapolation from the control region to the signal region. For 
the WW background in the A^jet < 1 categories and the top-quark background in the A^jet - 1 and A^jet > 2 
jet categories, theoretical uncertainties on the extrapolation are evaluated according to the prescription of 
Ref. [88]. The uncertainties include the impact of missing higher-order QCD corrections, PDF variations 
and MC modelling. For backgrounds normalised to the theoretical prediction without use of a control 
region, a similar prescription is followed for the acceptance uncertainty, and uncertainties on the predicted 
inclusive cross section also apply. 

For the WW background, the uncertainties on the control region extrapolation in the 0- and 1-jet categories 
amount to 4.2% and 9.7%, respectively. In the > 2 jet category, the WW yield is taken from the theoretical 
expectation. The PDF uncertainty on the cross sections of the qq —> WW + 2 jets and gg WW + 2 jets 
processes is evaluated to be 4% and that on the acceptance of these processes to be 2%. The QCD scale 
uncertainty on QCD WW + 2 jets is 14%, while that on the acceptance is 20%. The modelling uncertainty 
is derived by comparing samples generated with Sherpa and MadCraph [91] generators. This uncertainty 
amounts to 34% for QCD WW + 2 jets and 7% for EW WW + 2 jets after the selection on ma- In all jet 
categories, a correction is applied to the WW background to take higher-order EW corrections to the cross 
section into account. A conservative uncertainty of 100% is applied to this correction. 

Eor the top-quark background estimation in the 0-jet category, a 3.9% theoretical uncertainty is assigned 
from the use of simulated top-quark events to model the ratio of the signal region jet veto efficiency 
to the square of the efficiency of the additional-jet veto in the A-tagged control region in 

Eq. 3). An additional 4.5% theoretical uncertainty is used on the efficiency trest on the remaining selection 
that defines the Ajet = 0 signal region, which is also derived from simulated top-quark events. The most 
important component of these uncertainties is the variation among predictions by different MC generator 
and parton shower algorithms. Smaller uncertainties attributable to the QCD scale choice, PDE model, 
tt-Wt interference, and the single-top cross section are also included. Eor the top-quark background in the 
1- and > 2 jet categories, the uncertainties are respectively 5.7% and 9.8%. 

The main uncertainty on the VF-i-jets and multijet background predictions arise from the extrapolation 
factors relating anti-identified to identified leptons. Eor IT-i-jets, a modelling uncertainty on the simulation- 
based correction applied to the Z-i-jets extrapolation factor of about 20% and the statistical uncertainty 
from the Z-i-jets data used to measure the extrapolation factor contribute in roughly equal proportions. 
The main uncertainty on the multijet contribution arises from an uncertainty on the modelling of the cor¬ 
relation between the extrapolation factors for two anti-identified leptons in the same event. 

In the same-flavour channel in the > 2 jet category, the Drell-Yan background is non-negligible. The 
background is estimated from simulation in this final state as in the other final states, but a 15% theoretical 
uncertainty is assigned to the cross section of the process using the total relative theoretical uncertainty on 
the Z -I- 2 jets cross-section prediction in the high-mjj region, following Ref. [92]. 

In addition to the uncertainties on the normalisation of backgrounds, experimental and theoretical uncer¬ 
tainties on the mj shape model used in the fit are considered. Simulated data are used for the mj shape 
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for signal and all backgrounds except for lk+jets. Only shape variations which are statistically signific¬ 
ant compared to the statistical uncertainty from the simulation model of the shape are included. For the 
background model, uncertainties due to f^-tagging efficiency, lepton identification, trigger, and isolation 
efficiency scale factors fall into this category, as does the uncertainty on the extrapolation factors for the 
IF-I- jets and multijet background estimates. For the ggF signal, the experimental systematic uncertainties 
from are treated as mj shape uncertainties. The theoretical uncertainty from the interference weight¬ 
ing are also included as a systematic uncertainty on the mj shape for the ggF CPS signal samples. 


7.4 Uncertainties specific to the H WW (vqq analysis 

Uncertainties in the H WW —> (vqq analysis are analogous to those in the H WW —> (v(v analysis, 
but different uncertainties are prominent. The absence of exclusive jet binning means that there are no 
additional theoretical uncertainties on the signal acceptance. Experimental uncertainties, particularly those 
relating to the modelling of jet energies, are more important. 

In addition to the common experimental uncertainties, the H —> WW (vqq results have an uncertainty 
arising from the energy resolution of large-/? jets. This uncertainty is determined by matching recon¬ 
structed jets in simulated events with their associated particle-level jets, and computing the ratio of the 
reconstructed energy (mass) to their true values as determined after parton showering and hadronisation. 
Previous ATLAS studies of large-/? jet energy/mass resolution indicate that the resolution in simulation 
can vary by up to 20% [93]. Based on these studies, a systematic uncertainty is estimated by smearing the 
jet energies by a factor corresponding to an increase of 20% in their resolution in bins of pj and \q\. 

The dominant uncertainty on the background modelling in the H —> WW —> (vqq analysis is that on the 
shape of the miyjj spectrum for the IF-i-jets and top-quark backgrounds. This uncertainty is 50% of the 
difference befween fhe reweighfed and nominal miyjj disfribufions, as described in Secfion 6.4.1, and is 
faken as a sysfemafic uncerfainfy on fhe shape of fhe background disfribufion. 

The sysfemafic uncerfainfy on mulfijef producfion is defermined by comparing fhe nominal background 
esfimafe wifh an esfimafe derived using differenf lepfon selection criferia including a reversed isolafion 
requiremenf. Uncerfainfies are derived for differenf producfion mechanisms and decay channels as a 
funcfion of fhe Higgs boson mass, and vary from abouf 10% fo 100% of fhe mulfijef background. 


8 Signal and background predictions compared to data 

8.1 The H -* WW -* tvtv analysis 

In Table 7, fhe expecfed signal and background in fhe H —> WW —> (v(v analysis are summarized and 
compared fo fhe number of dafa evenfs passing fhe signal region selecfion. The CPS scenario is used for 
fhe signal, and fhe SM Higgs boson cross secfion is used fo normalise if. For fhis comparison only, fo 
give an indicafion of fhe sensifivify of fhe shape fif, an mj window is added fo fhe selecfion. The mj 
requiremenfs are chosen fo be abouf 80% efficienf for fhe signal, and are: ISOGeV < rriT < 270 GeV 
for niH - 300 GeV, 250 GeV < mj < 500 GeV for mh = 600 GeV, and 300 GeV < mj < 750 GeV for 


30 


Table 7 : Summary of the expected signal and background in the H —» WW Mv signal regions. The top table com¬ 
pares the observed number of candidate events in data A^data with the expected signal A^sig for several mn values and 
the total background A^bkg, along with its statistical (-y/A^bkg) and systematic (5A^bkP uncertainties. The predictions 
are quoted for an nij interval which is about 80% efficient for the signal (details in the text). The different-flavour 
and same-flavour final states are combined, and in the top table the A^jet = 0 and A^jet = 1 categories are summed. The 
bottom table shows the composition of the background for the analysis categories of different jet multiplicity. The 
VV background category includes all diboson processes except for WW. The Wj + jj column contains the sum of 
the ITH-jets and multijet backgrounds. 


Category 

CPS signal expectation 

niH 

^ Slg Slg 

Bkg. expectation 
fVbkg± 

Observed 

f^data 

fVjet < 1 

300 GeV 

144 

10 

961 

±31 

±47 

951 


600 GeV 

29 

3 

584 

±24 

± 12 

538 


900 GeV 

5 

2 

325 

± 18 

± 8 

290 

A^jet>2 

300 GeV 

3.1 

13.5 

18 

±4.2 

±5 

20 


600 GeV 

0.8 

4.2 

9.5 

± 3.1 

± 1.9 

15 


900 GeV 

0.2 

2.3 

5.6 

± 2.4 

± 1.5 

10 


Category 

ITIT 

VV 

top quark 

Wj + jj 

Z/y* 

H[125 GeV] 

fVjet = 0 

54.8% 

3.5% 

37.3% 

2.6% 

1.6% 

0.1% 

fVjet=l 

40.7% 

3.5% 

49.3% 

3.4% 

2.9% 

0.2% 

Vjet > 2 

24.6% 

2.3% 

36.3% 

2.5% 

30.7% 

3.6% 


mu - 900 GeV. The predicted event yields given do not include adjustments resulting from the fit to 
data. 

The systematic uncertainty on the background is derived from the expected uncertainty on the fit results, 
by fitting the nominal expected signal and background using the profile likelihood described in Section 9. 
The total expected uncertainty on the fit result has contributions from the theoretical uncertainties on the 
signal and background predictions, the experimental uncertainties, and the statistical uncertainties. After 
removing the signal theoretical uncertainties and the statistical uncertainty in quadrature, the approximate 
systematic uncertainty on the background, accounting for all systematic correlations among the back¬ 
grounds, can be extracted. 

The fractional background composition for each Ajet category is given in the bottom half of the table. 
The relative contributions of the various sources to the total background does not vary substantially as a 
function of nix, so the mx selection described above is not applied. 

Figure 5 shows mx distributions in the signal regions, separately for the different- and same-flavour chan¬ 
nels and for each jet category. No significant data excess is observed in any final state. The deficit of data 
events in the high-nix region of the 0-jet different-flavour final state has been investigated and no underly¬ 
ing systematic experimental or modelling effect has been identified. In particular, no correlated kinematic 
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effects are observed in the data, an alternative WW MC event generator (MC@NLO [94]) does not qual¬ 
itatively change the level of disagreement, and the mj distribution in the WW control region (Figure 1) 
does not show a comparable deficit. 


8.2 The H WW (vqq analysis 

In Table 8, the expected signal and background in the H WW tvqq analysis are summarized and 
compared to the number of data events passing the signal region selection. For this comparison only, 
to give an indication of the sensitivity of the shape fit, a requirement that mivjj be close to the mass 
hypothesis is added to the m//-dependent signal region selection. The requirement is \mtyjj - m//| < 
200 GeV), except for mn = 300 GeV, for which the lower bound is 200 GeV, corresponding to the lower 
ni{yjj considered in the analysis. The systematic uncertainty on the background is derived by fitting the 
nominal expected signal and background, following the procedure used for the H WW —> analysis. 

The fractional background composition for the ggF and VBF analysis categories and the ^//-dependent 
selection is shown in the bottom half of the table. The relative contributions of the various sources to the 
total background do not vary substantially as a function of m^yjj, so the ifiiyjj window described above is 
not applied. 

Figure 6 show the distributions and the ratio of data to background expectation for the WCR (top), 
TopCR (middle), and SR (bottom) after the ggF preselection on the left hand side of the figure. Shown on 
the right are the corresponding distributions for the VBF preselection. These distributions do not include 
the background normalisations applied by the fit to the control regions, but the miyjj reweighting described 
in Section 6.4.1 is applied. 

The Higgs boson signal yield in each final state is determined using a binned maximum likelihood fit 
to the observed niiyjj distribution in the range 200 GeV < m^yjj < 2000 GeV. For the mn - 500 GeV 
selection, the control and signal regions distributions are shown in Figure 7. In these distributions, the 
miyjj reweighting is applied and the background normalisations are corrected using the results of the fit to 
the signal and control regions. There is no indication of a significant excess of data above the background 
expectation. A slight deficit can be seen in the VBF channel in the centre and lower panels of Figure 6 but 
its effect is mitigated by the mass-dependent selection. 


9 Results and Interpretations 

9.1 Statistical methodology 

The methodology used to derive statistical results is described in detail in Ref. [95]. A likelihood function 
X is defined using the distributions of the discriminant for events in the signal region of each analysis 
category, namely, the 0-, 1- and > 2-jet categories in the H WW —» Mv analysis and the ggF and VBF 
categories in the H WW —> ivqq analysis. The likelihood is a product of Poisson functions over the 
bins of the discriminant in the signal regions and ones describing the total yield in each control region. 
Each systematic uncertainty is parameterised by a corresponding nuisance parameter 6 modelled by a 
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Figure 5; Transverse mass distributions in the 0-jet (top), 1-jet (middle) and > 2 jet (bottom) categories in the 
H —> WW —> signal regions. Different-flavour (left) and same-flavour (right) final states are shown. The 

histograms for the background processes, including the observed H[125 GeV], are shown stacked, and the distribu¬ 
tion for a hypothetical CPS signal process with nin - 600 GeV and the cross section predicted by the SM for that 
mass is overlaid. The combined statistical and systematiJTincertainties on the prediction are shown by the hatched 
band in the upper pane and the shaded band in the lower pane. In each figure, the last bin contains the overflow. 
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Figure 6; Distributions of the invariant mass m(yjj of the WW system and the ratio of data to background expectation 
for the WCR (top), TopCR (middle), and SR (bottom) after the ggF preselection (left) and the VBF preselection 
(right) for the H WW —> ivqq analysis. The histograms for the background processes are shown stacked, and the 
distribution for a hypothetical CPS signal process with nin - 500 GeV and the cross section predicted by the SM for 
that mass is overlaid. All the flavour and charge categories are summed in each plot. No normalisation scale factors 
are applied to the top-quark or W background samples. The combined statistical and systematic uncertainties on the 
prediction are shown by the hatched band in the upper pane and the shaded band in the lower pane. In each figure, 
the last bin contains the overflow. 
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Figure 7; Distributions of the invariant mass m(yjj of the WW system and the ratio of data to background expectation 
for the WCR (top), TopCR (middle), and SR (bottom) after the ggF 500 GeV selection (left) and the VBF 500 GeV 
selection (right) in the H —> WW (vqq analysis. These plots are after fitting for the nin - 500 GeV hypothesis, 
and the binning used is identical to the binning used in the fit. The histograms for the background processes are 
shown stacked, and the distribution for a hypothetical Cl^ signal process with nin - 500 GeV and the cross section 
predicted by the SM for that mass is overlaid. All the flavour and charge categories are summed in each plot. The 
combined statistical and systematic uncertainties on the prediction are shown by the hatched band in the upper pane 
and the shaded band in the lower pane. In each figure, the last bin contains the overflow. 

























































































Table 8: Summary of the expected signal and background in the H —> WW —> tvqq signal regions. The top table 
compares the observed number of candidate events in data A^data with the expected total signal, Nsig, for several nin 
values and the total background A^bkg, along with its statistical (-y/A^bkg) and systematic (AA^bkP uncertainties. The 
predictions are quoted for an mcyjj interval around the nin hypothesis (details in the text). The bottom table shows 
the composition of the background for the ggF and VBF analysis categories at several nin hypotheses. The “other” 
column contains the Z+jets and diboson backgrounds. 



CPS signal expeclalion 

Bkg. expeclalion 

Observed 

Category 

mu 

Slg 

A^VBF 

Slg 

f^bkg 

± V-^bkg ± dVbkg 

f^data 

ggF 

300 GeV 

1320 

100 

112890 

±340 ±460 

111199 


600 GeV 

440 

40 

23680 

± 150 ± 131 

23397 


900 GeV 

40 

10 

1940 

± 40 ± 56 

1754 

VBF 

300 GeV 

24 

41 

2282 

± 48 ± 33 

2090 


600 GeV 

24 

34 

850 

± 29 ± 19 

829 


900 GeV 

3 

11 

52 

± 7 ± 7 

68 


Category 

mu 

VF±jels 

lop 

mullijel 

olher 

ggF 

300 GeV 

70% 

20% 

3% 

7% 


600 GeV 

71% 

19% 

2% 

7% 


900 GeV 

73% 

14% 

3% 

10% 

VBF 

300 GeV 

58% 

33% 

5% 

4% 


600 GeV 

61% 

27% 

6% 

5% 


900 GeV 

52% 

32% 

4% 

12% 


Gaussian function (the set of all such nuisance parameters is 6). The modification of affected event yields 
is implemented as a log-normal distribution parameterised by 9 to prevent predicted event yields from 
taking unphysical values. 

In the H WW iv€v analysis, the htt distributions in the signal regions are divided into ten, six and 
four bins, respectively, for A^jet = 0, 1 and > 2. The bins are of variable widths such as to have the same 
number of expected signal events in each bin. Because the signal is peaked in mx and the background is 
not, this binning strategy improves the sensitivity by producing bins with different signal-to-background 
ratios and is robust against statistical fluctuations in the background model. 

In the H WW £vqq final slate, Ihe slalislical analysis is performed using a variable hi range and 
number of bins, adapted fo each mu hypofhesis and production mode. The bins are always of equal 
widfh, because fhe me^jj disfribufion for signal peaks more sfrongly fhan fhe signal mj disfribufion in 
fhe H WW iv€v channel, so fhe discriminanf produces bins wifh sufiicienlly differenl signal-fo- 
background ralios wifhouf furlher opfimisafion. The bin widfh is chosen lo ensure adequafe slalislical 
precision for fhe background predictions. The search is always preformed in a 700 GeV-wide window in 
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mivjj, enclosing the resonance peak and as much of the tails as feasible for non-NWA signal models. The 
range 200-900 GeV is used for the 300 GeV hypothesis, and 500-1200 GeV for the 1000 GeV hypothesis. 
The ggF (VBF) search uses 35 (17) bins for the 300 GeV hypothesis. The number of bins decreases with 
increasing niH, and 12 (6) bins are used for the ggF (VBF) search for niH - 1000 GeV. Mass hypotheses 
above 1 TeV are also tested, but the binning and fit range for the 1 TeV mass hypothesis are maintained, 
because there is insufficient data and simulated events to populate the background model at higher values 
of rriiyjj. 

Both the H WW Mv and H WW —> (vqq analyses have signal regions optimised for the VBF 
and ggF signal production modes, but the presence of both signal processes is accounted for in all signal 
regions. Limits are obtained separately for ggF and VBF production in all interpretations. To derive the 
expected limits on the ggF production mode, the VBF production cross section is set to zero, so that the 
expected limits correspond to the background-only hypothesis. To derive the observed limits on the ggF 
(VBF) production mode, the VBF (ggF) production cross section is treated as a nuisance parameter in 
the fit and profiled using a flaf prior, as is used for fhe normalisafion of backgrounds using CRs. This 
approach avoids making any assumpfion on fhe presence or absence of fhe signal in ofher producfion 
modes, by using fhe signal regions fhemselves fo set the normalisation of the production mode not being 
tested. 

The modified frequenfisf mefhod known as CLg, combined wifh fhe asympfofic approximafion, is used fo 
compufe 95% CL upper limifs [96, 97]. The mefhod uses a fesf sfafisfic a funcfion of fhe signal sfrengfh 
yU which is defined as fhe rafio of fhe measured cru x BR(// ^ WW) fo fhaf predicfed.^ The fesf sfafisfic 
is defined as: 

qf,^-l\n{£{ii-\)l£{fi-e)) (9) 

The denominator does nof depend on /r. The quanfifies jl and 0 are fhe values of yu and 6, respectively, fhaf 
unconditionally maximise £. The numerafor depends on fhe values Ofj fhaf maximise £ for a given value 
of yU. 


9.2 Upper limits from the H WW (v(v analysis 

Figure 8 shows fhe 95% CL upper limifs on cth x BR(// ^ WW) as a funcfion of mh for a Higgs boson in 
fhe CPS scenario, separafely for ggF and VBF producfion, in fhe mass range 220 GeV < mn < 1000 GeV. 
Figure 9 shows fhe upper limifs on a Higgs boson wifh a narrow widfh in fhe range 300 GeV < mh < 
1500 GeV, separately for ggF and VBF producfion. Below 300 GeV, fhe limifs in fhe NWA scenario 
are expecfed fo be similar fo fhe CPS scenario, as fhe widfh is small enough in fhe laffer case fo have a 
negligible effecf. 

The sysfemafic uncerfainfies wifh fhe largesf effecf on fhe observed limifs af mu = 300 GeV, in approx- 
imafe order of imporfance, are fhose related fo fhe modelling of fhe WW background, fhe f^-jef fagging 
efficiency, fhe jef energy scale and resolufion, fhe fop-quark background modelling in fhe Vjet = 0 caf- 
egory, fhe QCD scale uncerfainfies on fhe signal from fhe exclusive jet multiplicity categories, and the jet 
energy scale and resolution. As the mass hypothesis increases, the experimental systematic uncertainties 


^ The SM cross-section prediction is used to define fi for the NWA and intermediate-width scenarios. 
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Figure 8; 95% CL upper limits on cth x BR(// WW) from the H WW —> Mv analysis for the CPS scenario. 
Limits for ggF production (left) and VBF production (right) are shown. The green and yellow bands show the +lcr 
and +2cr uncertainties on the expected limit. 


diminish in relative importance and the leading sources of uncertainty are on the signal cross section in 
exclusive A^jet categories and the WW background model. Also, for the CPS scenario, the uncertainty on 
the interference weighting becomes important at high mn- 

The relative importance of systematic and statistical uncertainties for the analysis can be illustrated by 
recalculating the limits with the systematic uncertainties omitted. For the CPS scenario, the ggF limits 
decrease by about 20% for mh = 300 GeV and by about 10% for mn = 1000 GeV. For high mass 
hypotheses, higher values of niT, where there are fewer events, are implicitly tested. Similarly, there are 
fewer candidates in the Ajet > 2 category, so that the limits are less sensitive to systematic uncertainties 
than the corresponding ggF limits. For the CPS VBF limit, removing the systematic uncertainties has 
about an 8% effect on the observed limit at mh = 300 GeV and a negligible effect at mh = 1000 GeV. 

The deficit at mx > 450 GeV observed in Figure 5 results in a stronger limit than that predicted for 
background-only in the ggF production mode in both signal models, although they are consistent within 
the given uncertainties. For signal, the relation nix < nin holds, so the observed limits are stronger than 
the expected ones above this threshold of about 450 GeV. 


9.3 Upper limits from the H WW (vqq analysis 

Limits are derived following the same procedure as for the H —> WW —> £v£v channel. Figure 10 shows 
the 95% CL upper limits on cth x BR(// —> WW) as a function of mu for the CPS scenario, separately 
for ggF and VBF production, in the mass range 300 GeV < mu < 1000 GeV. The limits derived from the 
H WW —> ivqq analysis are comparable to those derived from the H —> WW —> Mv analysis. Figure 11 
shows the upper limits on a Higgs boson with a narrow width in the range 300 GeV < mu < 1500 GeV, 
separately for ggF and VBF production. 
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Figure 9; 95% CL upper limits on cth x BR(H —> WW) from the H —> WW Mv analysis for a signal with a 
narrow width. Limits for ggF production (left) and VBF production (right) are shown. The green and yellow bands 
show the +lcr and ±2cr uncertainties on the expected limit. 


The systematic uncertainties with the largest effect on the observed limits at mn = 300 GeV are those 
related to the jet energy scale and resolution, the multijet background estimation, and the ft-jet tagging, 
particularly the uncertainty on the rate for mistakenly tagging a light-quark jet as a ^-jet. As the tested 
mass hypothesis increases, the mtyjj shape uncertainties on the VF-t jets and top backgrounds become the 
leading sources of uncertainty and the multijet background systematic uncertainties become negligible. 
For the CPS scenario, the uncertainty on the interference weighting also becomes a dominant systematic 
uncertainty at high mn- 

To show the overall effect of systematic uncertainties, the exercise done for the H WW —> Mv analysis 
is repeated. If systematic uncertainties are omitted, the observed ggF limits in the CPS scenario decrease 
by 66% at niH = 300 GeV and 40% at mh - 1000 GeV. The corresponding VBF limits decrease by about 
40% and 20%, respectively. The trends relative to mh and ggF vs. VBF are similar to what is seen in the 
dilepton final state, but systematic uncertainties have a larger effect on the H WW £vqq limits than 
on the H WW —> limits because the larger candidate event samples in the former analysis result in 

smaller statistical uncertainties. 

The downward excursions of the observed limits compared to the expected ones seen for mu > 600 GeV in 
Figure 10 in the ggF category and for tue » 750 GeV in both categories in Figure 11 have been investigated 
and no underlying systematic effect identified. In particular, the simultaneous dip in the ggF and VBF 
NWA limits at mn ~ 750 GeV is attributable to a coincidence of deficits in the data in the statistically 
independent ggF and VBF SRs at that value of mtyjj. 


9.4 Combined upper limits 

This section presents 95% CL upper limits on the production of high-mass Higgs bosons in the CPS and 
NWA scenarios from a combination of the H WW —> Mv and H —> WW Cvqq final states. In 
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Figure 10; 95% CL upper limits on cth x BR(H WW) from the H —> WW —> fvqq analysis for the CPS scenario. 
Limits for ggF production (left) and VBF production (right) are shown. The green and yellow bands show the +lcr 
and +2cr uncertainties on the expected limit. 




Figure 11: 95% CL upper limits on cth x BR(// —> WW) from the H —> WW —> fvqq analysis for a signal with a 
narrow width. Limits for ggF production (left) and VBF production (right) are shown. The green and yellow bands 
show the +lcr and +2cr uncertainties on the expected limit. 


the statistical combination, the likelihood function is constructed from the signal and background prob¬ 
ability density functions from the two analyses. The combination takes into account all statistical and 
systematic uncertainties in both analyses. In particular, correlated effects of given sources of systematic 
uncertainties in the two final slates are taken into account correctly. These correlated effects arise from 
sources of uncertainty common to the final sfales, for example, those related to detector response affecting 
the reconstruction, identification and calibration of electrons, muons, jets, and Z>-tagging, as well 
as the integrated luminosity. Systematic uncertainties that affect both final stales are correlated in the 
combination unless there is a specific reason nof fo correlale Ihem. 
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Figure 12: 95% CL upper limits on cth x BR(// —> WW) for the CPS scenario from the combination of the 
H —> WW Ivtv and H —> WW —> tvqq final states. Limits for ggF production (left) and VBF production 
(right) are shown. The green and yellow bands show the +lcr and +2cr uncertainties on the expected limit. 


Since the H WW —> £vqq analysis sets upper limits starting at a Higgs boson mass hypothesis of 300 
GeV, the combination is performed starting at rtiH = 300 GeV. In the CPS scenario, the upper range of the 
combination is mn = 1000 GeV since neither analysis performs the search above this mass because of the 
large width. In the NWA case, the upper range of the combination extends to mh = 1500 GeV. 

Figure 12 shows combined upper limits separately on the ggF and VBF production modes for a Higgs 
boson in the CPS scenario. As in the case of the H —> WW —> Mv and H —> WW —> ivqq final states, 
when expected limits on a given production mode are extracted, the cross section of the other production 
mode is set to zero, while for deriving observed limits on the production mode, the cross section of the 
other mode is profiled using dafa. 

Figure 13 shows fhe limits on cth x BR(// ^ WW) as a function of mh for a narrow-width Higgs boson, 
separately for the ggF and VBF production modes. As in the CPS scenario, when observed limits on a 
given production mode are extracted, the strength parameter of the other production mode is profiled as a 
nuisance paramefer in the fit. 


9.5 Results in the intermediate-width scenario 

The data can also be interpreted in terms of an additional Higgs boson with a width intermediate between 
the narrow-width approximation and the CPS scenario. This interpretation is motivated by the electroweak 
singlet (EWS) model, and assumes that the production cross sections and partial widths of the heavy Higgs 
boson are related to those of the SM Higgs boson by a single, constant scale factor This allows 

combination of the ggF and VBF production modes as well as the two WW decay channels considered 
here. Non-SM decay modes, possible in the EWS model, are not considered in this analysis and the 
branching ratios of the heavy Higgs boson are the same as for a hypothetical SM Higgs boson of the same 


41 











Figure 13: 95% CL upper limits on cth x BR(i/ WW) for a signal with a narrow width from the combination 
of the H —» WW —» Mv and H —> WW (vqq final states. Limits for ggF production (left) and VBF production 
(right) are shown. The green and yellow bands show the +lcr and +2cr uncertainties on the expected limit. 


mass. The cross section, width and branching ratio of the heavy Higgs boson can be expressed as follows: 

(TH = X<Th,sm 

Th = X Fz/ sm (10) 

BR, = BRsm.i- 

where cr//,sM> Fand BRsm,/ are the cross section, total width and branching ratio to decay mode i 
of a SM Higgs boson with mass mn, respectively. The parameters of a true electroweak singlet model 
are substantially constrained by measurements of the Higgs boson at mn » 125 GeV [3]. The treatment 
described here allows a greater spectrum of possible widths to be explored. 

Figure 14 shows upper limits in the H WW —> £v£v channel as a function of mu in the intermediate- 
width scenario for widths in the range 0.2F//,sm < F// < 0.8F//sm- Limits are shown on cth x BR(// ^ 
WW) divided by k'^ to facilitate readability, since otherwise the limit curves corresponding to the vari¬ 
ous k'^ values approximately coincide. This feature indicates that the H —> WW £v(v channel has 
little sensitivity to the width of the resonance. Similarly, Figure 15 shows corresponding limits for 
the H WW —> ivqq analysis, and Figure 16 shows the combination of the H WW £v£v and 
H WW £vqq analyses. 
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Figure 14: 95% CL upper limits in the intermediate-width scenario on ^ x BR(// —» WW) from the 
H —» WW —> Mv analysis for a heavy scalar resonance with width in the range 0.2rH,sM < F// < 0.8Fh,sm- 
The ggF and VBF production modes have been combined. 



Figure 15: 95% CL upper limits in the intermediate-width scenario on ^ x BR(// —» WW) from the 
H —» WW —» tvqq analysis for a heavy scalar resonance with width in the range 0.2F//,sm < Fh < 0.8Fh,sm- 
The ggF and VBF production modes have been combined. 
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Figure 16: 95% CL upper limits in the intermediate-width scenario from the combination of the H —> WW —» Mv 
and H —» WW —> ivqq hnal states. Limits are set on ^ x BR(H —» WW) for a heavy scalar resonance with a width 
in the range 0.2r/f,sM < F// < 0.8 F//,sm- The ggF and VBF production modes have been combined. 


10 Conclusion 

A search for a high-mass Higgs boson in the H —> WW —> Mv and H WW —> ivqq channels by the 
ATLAS experiment at the LHC is reported. The search uses proton-proton collision data at a centre-of- 
mass energy of 8 TeV corresponding to an integrated luminosity of 20.3 fb“^ No excess of events beyond 
the Standard Model background prediction is found. Upper limits are set on the product of the production 
cross section and the H WW branching ratio in three different scenarios: a high-mass Higgs boson 
with a CPS lineshape and the width predicted for a SM Higgs boson, one with a narrow width, and one 
with an intermediate width. 

For all mass hypotheses tested, the strongest upper limits are observed for the narrow-width scenario. At 
mu = 300 GeV, these are cru x BR{H WW) < 830 fb at 95% CL for the gluon-fusion production 
mode and cru x BR(// —> WW) < 240 fb at 95% CL for the vector-boson fusion production mode. For 
mu = 1500 GeV, the highest mass-point tested, the cross section times branching ratio is constrained to 
be less than 22 fb and 6.6 fb at the 95% CL for the ggF and VBF production modes, respectively. 

The limits in the CPS and intermediate-width scenarios are qualitatively similar but somewhat weaker 
due to the increased resonance width, particularly for mu approaching 1 TeV. For the CPS scenario, the 
combined 95% CL upper limits on cru x BR(// ^ WW) for the ggF production mode range from 990 fb 
at mu - 300 GeV to 35 fb at mu = 1000 GeV. For the VBF production mode, the equivalent values are 
230 fb and 27 fb, respectively. These results are a substantial improvement over the previous results from 
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the ATLAS experiment in terms of both the range of mu explored and the cross section times branching 
ratio values excluded. 
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